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Abstract 

We compute the chromo-field distributions of static color dipoles in the 
fundamental and adjoint representation of SU{Nc) in the loop-loop correla¬ 
tion model and find Casimir scaling in agreement with recent lattice resnlts. 
Our model combines perturbative gluon exchange with the non-perturbative 
stochastic vacuum model which leads to confinement of the color charges in 
the dipole via a string of color fields. We compnte the energy stored in the con¬ 
fining string and use low-energy theorems to show consistency with the static 
quark-antiquark potential. We generalize Meggiolaro’s analytic continuation 
from parton-parton to gange-invariant dipole-dipole scattering and obtain a 
Euclidean approach to high-energy scattering that allows us in principle to 
calculate 5-matrix elements directly in lattice simulations of QCD. We apply 
this approach and compute the 5-matrix element for high-energy dipole-dipole 
scattering with the presented Euclidean loop-loop correlation model. The re¬ 
sult confirms the analytic continuation of the gluon field strength correlator 
used in all earlier applications of the stochastic vacunm model to high-energy 
scattering. 
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1 Introduction 


The structure of the QCD vacuum is responsible for color conhnement, spontaneous 
chiral symmetry breaking, and dynamical mass generation [1], Hadronic reactions 
are expected to show further manifestations of a non-trivial QCD vacuum. It is 
indeed a key issue to unravel the effects of conhnement and topologically non-trivial 
gauge held conhgurations (such as instantons) on such reactions [2,3]. Moreover, it 
would be a signihcant breakthrough to understand the size, behavior and growth of 
hadronic cross sections with increasing c.m. energy from the QCD Lagrangian. 

Lattice QCD is the principal theoretical tool to study the QCD vacuum from 
hrst principles. Numerical simulations of QCD on Euclidean lattices give strong 
evidence for color conhnement and spontaneous chiral symmetry breaking and de¬ 
scribe dynamical mass generation from the QCD Lagrangian [4-6]. However, since 
lattice QCD is limited to the Euclidean formulation of QCD, it cannot be applied in 
Minkowski space-time to simulate high-energy reactions in which particles are mov¬ 
ing near the light cone. Furthermore, although lattice investigations have signih- 
cantly enhanced our understanding of non-perturbative phenomena and particularly 
conhnement, one can quote the concluding sentence of Greensite’s recent review [7]: 
“The conhnement problem is still open, and remains a major intellectual challenge 
in our held.” Here (phenomenological) models that allow analytic calculations are 
important as they provide valuable complementary insights. 

In this work we introduce the Euclidean version of the loop-loop correlation model 
(LLCM) which has been developed in Minkowski space-time to describe high-energy 
reactions of hadrons and photons [8] on the basis of a functional integral approach [9- 
13]. The central element in our approach is the gauge-invariant Wegner-Wilson 
loop [14,15]: The physical quantities considered are obtained from the vacuum ex¬ 
pectation value (VEV) of one Wegner-Wilson loop, {Wr[C]), and the correlation of 
two Wegner-Wilson loops, [Ci]H72 [C 2 ]). Here r(j) indicates the SUi^NQ rep¬ 
resentation of the Wegner-Wilson loops which we keep as general as possible. We 
express {Wr[C]) and (H 7 i[Ci]H 72 [^ 2 ]) in terms of the gauge-invariant bilocal gluon 
held strength correlator integrated over minimal surfaces by using the non-Abelian 
Stokes theorem [16] and a matrix cumulant expansion [17] in the Gaussian approx¬ 
imation. The latter approximation relies on the assumption of a Gaussian domi¬ 
nance in the correlations of gauge-invariant non-local gluon held strengths, i.e. the 
dominance of the bilocal correlator over higher ones, and is supported by lattice 
investigations [18]. In our model this Gaussian approximation leads directly to the 
Casimir scaling of the static quark-antiquark potential which for SU (3) has clearly 
been conhrmed on the lattice [19,20]. We decompose the gauge-invariant bilocal 
gluon held strength correlator into a perturbative and a non-perturbative compo- 
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nent: The stochastic vacuum model (SVM) [21] is used for the non-perturbative 
low-frequency background held and perturbative gluon exehange for the additional 
high-frequency contributions. This combination allows us to describe long and short 
distance correlations in agreement with lattice calculations of the gluon held strength 
correlator [18,22-24]. Moreover, it leads to a static quark-antiquark potential with 
color Coulomb behavior for small and conhning linear rise for large source separa¬ 
tions. We calculate the static quark-antiquark potential with the LLCM parameters 
determined in hts to high-energy scattering data [8] and hnd good agreement with 
lattice data. We thus have one model that describes both static hadronic prop¬ 
erties and high-energy reactions of hadrons and photons in good agreement with 
experimental and lattice QCD data. 

We apply the LLCM to compute the chromo-electric fields generated by a static 
color dipole in the fundamental and adjoint representation of SU{Nc). The non- 
perturbative SVM component describes the formation of a color hux tube that con- 
hnes the two color sources in the dipole [25] while the perturbative component leads 
to color Coulomb helds. We hnd Casimir scaling for both the perturbative and non- 
perturbative contributions to the chromo-electric helds again as a direct consequence 
of the Gaussian appoximation in the gluon held strengths. The mean squared radius 
of the conhning QCD string is calculated as a function of the dipole size. Transverse 
and longitudinal energy density prohles are provided to study the interplay between 
perturbative and non-perturbative physics for diherent dipole sizes. The transition 
from perturbative to string behavior is found at source separations of about 0.5 fm 
in agreement with the recent results of Liischer and Weisz [26]. 

The low-energy theorems, known in lattice QCD as Michael sum rules [27], relate 
the energy and action stored in the chromo-helds of a static color dipole to the 
corresponding ground state energy. The Michael sum rules, however, are incomplete 
in their original form [27]. We present the complete energy and action sum rules [28- 
30] in continuum theory taking into account the contributions to the action sum rule 
found in [31] and the trace anomaly contribution to the energy sum rule [28]. Using 
these low-energy theorems, we compare the energy and action stored in the conhning 
string with the conhning part of the static quark-antiquark potential. This allows 
us to conhrm consistency of the model results and to determine the values of the 
Callan-Symanzik function and the strong coupling at the renormalization scale 
at which the non-perturbative SVM component is working. The values obtained for 
and as are compared to model independent QCD results for the Callan-Symanzik 
function. Earlier investigations along these lines have been incomplete since only 
the contribution from the traceless part of the energy-momentum tensor has been 
considered in the energy sum rule. 

To study the ehect of the conhning QCD string examined in Euclidean space- 
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time on high-energy reactions in Minkowski space-time, an analytic continuation 
from Euclidean to Minkowski space-time is needed. For investigations of high-energy 
reactions in our Euclidean model, the gauge-invariant bilocal gluon held strength 
correlator can be analytically continued from Euclidean to Minkowski space-time. 
This analytic continuation has been introduced for applications of the SVM to 
high-energy reactions [11-13] and is used in our Minkowskian applications of the 
LLCM [8,32-34]. Recently, an alternative analytic continuation for parton-parton 
scattering has been established in the perturbative context by Meggiolaro [35]. This 
analytic continuation has already been used to access high-energy scattering from 
the supergravity side of the AdS/CFT correspondence [36], which requires a posi¬ 
tive dehnite metric in the dehnition of the minimal surface [37], and to examine the 
effect of instantons on high-energy scattering [38]. 

In this work we generalize Meggiolaro’s analytic continuation [35] from parton- 
parton to gauge-invariant dipole-dipole seattering such that S'-matrix elements for 
high-energy reactions can be computed from conhgurations of Wegner-Wilson loops 
in Euclidean space-time and with Euclidean functional integrals. This shows how 
one can access high-energy reactions directly in lattice QCD. First attempts in this 
direction have already been carried out but only very few signals could be extracted, 
while most of the data was dominated by noise [39]. We apply this approach to 
compute the scattering of dipoles at high-energy in the Euclidean LLCM. We recover 
exactly the result derived with the analytic continuation of the gluon held strength 
correlator [8]. This conhrms the analytic continuation used in all earlier applications 
of the stochastic vacuum model to high-energy scattering [11-13,40-48] including 
the Minkowskian applications of the LLCM [8,32-34]. In fact, the S'-matrix element 
obtained has already been used as the basis for a unihed description of hadronic 
high-energy reactions [8], to study saturation effects in hadronic cross sections [8, 
32,34], and to investigate manifestations of the conhning QCD string in high-energy 
reactions of photons and hadrons [33]. 

The outline of the paper is as follows. In Sec. 2 the LLCM is introduced in its 
Euclidean version and the general computations of {Wr[C]) and (Wn[C'i]lTr. 2 [C 2 ]) 
are presented. Based on these evaluations, we compute the potential of a static 
color dipole in Sec. 3 and the associated chromo-held distributions in Sec. 4 with 
emphasis on Casimir scaling and the interplay between perturbative color Coulomb 
behavior and non-perturbative formation of the conhning QCD string. In Sec. 5 low- 
energy theorems are discussed and used to show consistency of the model results 
and to determine the values of (3 and at the renormalization scale at which the 
non-perturbative SVM component is working. In Sec. 6 the Euclidean approach to 
high-energy scattering is presented and applied to compute high-energy dipole-dipole 
scattering in our Euclidean model. In the Appendixes we review the derivation of 
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the non-Abelian Stokes theorem, give parametrizations of the loops and the minimal 
surfaces, and provide the detailed computations for the results in the main text. 


2 The Loop-Loop Correlation Model 

In this section the vacuum expectation value of one Wegner-Wilson loop and the 
correlation of two Wegner-Wilson loops are computed for arbitrary loop geometries 
within a Gaussian approximation in the gluon held strengths. The results are applied 
in the following sections. We describe our model for the QCD vacuum in which 
the stochastic vacuum model [21] is used for the non-perturbative low-frequency 
background held (long-distance correlations) and perturbative gluon exchange for 
the additional high-frequency contributions (short-distance correlations). 

2.1 Vacuum Expectation Value of one Wegner-Wilson Loop 

A crucial quantity in gauge theories is the Wegner-Wilson loop operator [14,15] 

IW[C'] = fr.Pexp dZ^^;](Z)f“ . (2.1) 

Concentrating on SU{Nc) Wegner-Wilson loops, where N^. is the number of colors, 
the subscript r indicates a representation of SU{Nc), Tr,. = Tr,,(- ■ ■ )/Trllr is the nor¬ 
malized trace in the corresponding color space with unit element 11^, S' is the strong 
coupling, and = ^“(Z)t“ represents the gluon held with the SU{Nc) group 

generators in the corresponding representation, that demand the path ordering 
indicated by V on the closed path C in space-time. A distinguishing theoretical fea¬ 
ture of the Wegner-Wilson loop is its invariance under local gauge transformations 
in color space. Therefore, it is the basic object in lattice gauge theories [4,14,15] 
and has been considered as the fundamental building block for a gauge theory in 
terms of gauge invariant variables [49]. Phenomenologically, the Wegner-Wilson 
loop represents the phase factor associated with the propagation of a very massive 
color source in the representation r of the gauge group SU{Nc). 

To compute the expectation value of the Wegner-Wilson loop (2.1) in the QCD 
vacuum 

(Wr[C]^^=(TTrVexp -tg £dZ,g;{Z)t<:: , (2.2) 

we transform the line integral over the loop C into an integral over the surface S 
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(2.3) 


with dS = C by applying the non-Abelian Stokes theorem [16] 

{Wr[C])^=(TTrVsexp -i^- Jda^,{Z)g;,{0,Z-Czo)K )^ , 

where Vs indicates surface ordering and O is an arbitrary reference point on the 
surface S. In Eq. (2.3) the gluon held strength tensor, Q^^^Z) = Q^^{Z)t^, is 
parallel transported to the reference point O along the path Czo, 

g,,{0, z- Czo) = MO, Z- Czo)GMZ)MO, Z; CzoM , (2.4) 

with the QCD Schwinger string 

MO,Z-,Czo)=Vexp -ig [ dZ^g;{Z)t: . (2.5) 

dCzo 

A more detailed explanation of the non-Abelian Stokes theorem and the associated 
surface ordering is given in Appendix A. 

The QCD vacuum expectation value {.. .)g represents functional integrals in 
which the functional integration over the fermion helds has already been carried out 
as indicated by the subscript G [10]. The model we use for the evaluation of (.. .)g is 
based on the quenched approximation which does not allow string breaking through 
dynamical quark-antiquark production. So far, it is not clear how to introduce 
dynamical quarks into this model. One suggestion is presented in Appendix A of 
Ref. [10]. 

Due to the linearity of the functional integral, (Tr^ • • •) = Trr{. • ■), we can write 

(Wr[C])^ = M(Vsexp jdaMZ)gMO,Z-,Czo)K )^ • ( 2 - 6 ) 

For the evaluation of (2.6), a matrix cumulant expansion'll used as explained in [10] 
(cf. also [17]) 

(Vs exp f da{Z)g{0,Z-,Czo) ) 

OO n 

= da(Xi)...da(XOifn(Xi,... ,X„)] , (2.7) 

n=l * ^ 

where space-time indices are suppressed to simplify notation. The cumulants Kn 
consist of expectation values of ordered products of the non-commuting matrices 
g{0, Z; Czo)- The leading matrix cumulants are 

Kgx) = {g{o,x-Cx))G, ( 2 . 8 ) 

R2(Xi,X2) = {Vs[g{0,XpCx,)g{0,X2,Cx,)])G 

-^[(^(0,Xi;CxJ)g(^(0,X2;CxJ)g + (1-2)] . (2.9) 
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Since the vacuum does not prefer a specific color direction, Ki vanishes and K 2 
becomes 

i^2(Xi,X2) = {Vs[Q{0,X^-Cx,)Q{0,X2,Cx,)])g ■ (2.10) 

Now, we approximate the functional integral associated with the expectation values 
(.. .)g as a Gaussian integral in the parallel transported gluon held strength (2.4). 
This Gaussian approximation is supported by lattice investigations [18] that show a 
dominance of the bilocal gauge-invariant gluon held strength correlator over higher- 
point non-local correlators. As a consequence of the Gaussian approximation, the 
cumulants factorize into two-point held correlators such that all higher cumulants, 
Kn with n > 2, vanish.^ Thus, (1 T).[G])g can be expressed in terms of K 2 

(M''r|Cl)^ = Tr^exp jda^,„{Xi) jda^^iX^) 

{Psl6;jO,Xr,Cx,o)K e‘,(0,A'2;Cx,o)itl)J (2.11) 

Due to the color neutrality of the vacuum, the gauge-invariant bilocal gluon held 
strength correlator contains a S function in color space, 

{£2 [6;jO,Xj-,Cx,o)6UO,X,_-,Cx,o)] =: A2,0;Cx.o,Cx,o) 

( 2 . 12 ) 

which makes the surface ordering Vs in (2.11) irrelevant. The tensor will be 

specihed in Sec. 2.3. With (2.12) and the quadratic Gasimir operator G 2 (r), 

V^V^=tl = C2ir)% , (2.13) 

Eq. (2.11) reads 

= fr^exp = exp , (2.14) 

where 

Xss ■= ^ Jdap^{Xi) Jdap^{X2)Fp^p^{Xi,X2,0-,CxiO,Cx2o)- (2-15) 

In this rather general result (2.14) obtained directly from the color neutrality of the 
QGD vacuum and the Gaussian approximation in the gluon held strengths, the more 
detailed aspects of the QGD vacuum and the geometry of the considered Wegner- 
Wilson loop are encoded in the function xss which is computed in Appendix G for 
a rectangular loop. 

^We are going to use the cumulant expansion in the Gaussian approximation also for perturba¬ 
tive gluon exchange. Here certainly the higher cumulants are non-zero. 
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In explicit computations we use for S the minimal surface, which is the planar 
surface spanned by the loop {C = dS) that leads most naturally to Wilson’s area 
law [21], Of course, the results should not depend on the choice of the surface. In 
our model the perturbative and non-perturbative non-conhning components satisfy 
this requirement. The non-perturbative conhning component in depends on 

the choice of the surface due to the Gaussian approximation and the associated 
truncation of the cumulant expansion. Since the minimal surface leads to a static 
quark-antiquark potential that is in good agreement with lattice data (see Sec. 3), 
we think that the minimal surface reduces the contribution from higher cumulants. 
Within bosonic string theory, our minimal surface represents the world-sheet of a 
rigid string: Our model does not describe fluctuations or excitations of the string 
and thus cannot reproduce the Liischer term which has recently been conhrmed by 
Liischer and Weisz [26]. 


2.2 The Loop-Loop Correlation Function 

The computation of the loop-loop correlation function (Wr][Ci\Wr 2 [C 2 \)G starts 
again with the application of the non-Abelian Stokes theorem [16] which allows us 
to transform the line integrals over the loops Gi 2 into integrals over surfaces S'! 2 
with dSi ^2 = Ci ^2 


WrACl]Wr,[C 2 ]) ={TTr,Vsexp 
G 


-^1 


X Tr^j exp 


^ J S2 


G 


(2.16) 


where Oi and O 2 are the reference points on the surfaces Si and S 2 , respectively, that 
enter through the non-Abelian Stokes theorem. In order to ensure gauge invariance 
in our model, the gluon held strengths associated with the loops must be compared 
at one reference point O. Due to this physical constraint, the surfaces Si and S2 
are required to touch at a common reference point Oi = O 2 = O. 

To treat the product of the two traces in (2.16), we transfer the approach of 
Berger and Nachtmann [45] (cf. also [ 8 ]) to Euclidean space-time. Accordingly, 
the product of the two traces, Tr,,j(-■ ■) Tr,, 2 (-■ ■), over SU^Nf) matrices in the 
ri and r 2 representations, respectively, is interpreted as one trace Trr^ 0 ,, 2 (-■ ■) ;= 
Tr,,^ 0 r 2 (' ■ ■ )/Trri(gir 2 (llri(gir 2 ) fh^t acts in the tensor product space built from the ri 
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and r 2 representations, 


G 


WrACl]Wr,[C2]) = {^Tr,^r,\\Vse^]^[-l^- j da^,{X,)g;,{0, X,; Cx,o) K^] 


11 , 


r2 


Si 


X 


11 


T1 


Vs exp [-*1 / dap^{X2)gl{0, X^; Cx,o) 4] 


'52 


G 


(2.17) 


With the identities 


exp(t“^) (g) 11^2 = exp(t“^ (g) llr-a) , (2-18) 

11^1 ® exp () = exp lXi^K^ ) , (2.19) 

the tensor products can be shifted into the exponents. Using the matrix multiplica¬ 
tion relations in the tensor product space 

® 11 , 2 ) ® 11 . 2 ) = 

® 11 , 2 ) ( 11.1 ® ty = 

and the vanishing of the commutator 

y^®ll,2,lln®4] =0i (2-21) 

the two exponentials in (2.17) commute and can be written as one exponential 


Vrit, ® 11^2 , 

y ® 4 1 (2-20) 


lU[C'i]lU[C'2 


G 


= (Tr 


r'i0r2 


P^exp 


. 9 


da^,{X)g^,{0,X;Cxo) 


G 


( 2 . 22 ) 


with the following gluon held strength tensor acting in the tensor product space: 


g^u{o,x-Cxo) 


^“,( 0 ,X;C'xo)(t“, ® 11,2) for X e , 
6 ;“,( 0 ,X;C'xo)(lln ® for X G ^2 • 


In Eq. (2.22) the surface integrals over Si and S 2 are written as one integral over 
the combined surface S' = -|- 5*2 so that the right-hand side (RHS) of (2.22) 

becomes very similar to the RHS of (2.3). This allows us to proceed analogously to 
the computation of {Wr[C])G in the previous section. After exploiting the linearity 
of the functional integral, the matrix cumulant expansion is applied, which holds 
for g^y{0,X-^Cxo) as well. Then, with the color neutrality of the vacuum and by 
imposing the Gaussian approximation now in the color components of the gluon held 
strength tensor,^ only the n = 2 term of the matrix cumulant expansion survives, 

^Note that the Gaussian approximation on the level of the color components of the gluon 
field strength tensor (component factorization) differs from the one on the level of the gluon 
field strength tensor (matrix factorization) used to compute {Wr\C]) in the original version of 
the SVM [21]. Nevertheless, with the additional ordering rule [25] explained in detail in Sec. 2.4 
of [50], a modified component factorization is obtained that leads to the same area law as the 
matrix factorization. 



which leads to 


WrACl]Wr,[C2[ 


= exp 


G 


(2.24) 


da^^iX,) dap^iX2)(Vs[g^A0, Xp Cx,o)GpaiO, X^; Cx,o)] 


G 


Using dehnition (2.23) and relations (2.20), we now redivide the exponent in (2.24) 
into integrals of the ordinary parallel transported gluon held strengths over the sep¬ 
arate surfaces Si and S 2 


fUri [Cl ] fUr2 [C2] ) = Tr^^ exp 

G 


(2.25) 


x<j y^da^.(Xi)y^da,.(X2)iP5[(^“.(0,Xi;Ux,o)^J.(0,X2;Cx2o))^(t“, ® <) 


>S 2 


+ da,AXi) da,^{X2) Vs fe,( 0 , Xp Ux,o)^J.( 0 , X2; Cx,o)) {K, ® 4 ) 


'S 2 


'Si 


G 


+ da,AXi) da,^{X2) Vs fc( 0 , Xp Ux,o)^J.( 0 , Xp Ux,o)) {K/n ® X2) 


'Si 


'Si 


G 


+ y^da^,(Xi)y^da,.(X2)iP5[(^“,(0,XpUx,o)^J.(0,X2;Ux,o))^(lln ® ^ 2 ^) 


Here the surface ordering Vs is again irrelevant due to the color neutrality of the 
vacuum (2.12), and (2.25) becomes 


HUi [Ci]Wr^ [C2] ) = Tr^i0r2 exp 

G 


XS1S2 + XS2S1 


(y ® o 


XsiSi 


{Kitri ® ^r-2) 


XS2S2 


(^n ® K 2 K 2 ) 


(2.26) 


with 


71 


XSiSj ^ I da,iy{Xi) / dapa{X2)F,i,,„{Xi,X2,0-,CxiO,Cx2o)- (2.27) 

V Si V Sj 


The symmetries in the tensor structure of [see Eqs. (2.42), (2.44), and (2.48)] 
lead to XS 1 S 2 = XS 2 Si- With the quadratic Casimir operator (2.13), our hnal Eu- 
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clidean result for general SU (iVc) representations ri and r 2 becomes^ 

(Wr,[Cl]Wr,[C2])^ (2.28) 

= frr,^r 2 exp - XS 1 S 2 {Ki ® ^“ 2 ) - XS 252 ) llri®r 2 

where Ilri 0 r 2 := ^ri ® 11^2. After specifying the representations ri and r 2 , the tensor 
product tri®r 2 := Ki ® K 2 ean be expressed as a sum of projection operators P* with 
the property P* ^^^ 0^2 = P* 

W, = E^'P- A. = A;fAhi^, (2,29) 

Tr^i(g)^2 (Pi j 

which corresponds to the decomposition of the tensor product space into irreducible 
representations. 


For two Wegner-Wilson loops in the fundamental representation of SU{Nc), ri = 
r 2 = Nc, which could describe the trajectories of two quark-antiquark pairs, the 
decomposition (2.29) becomes trivial: 


® - 2W * 2W “ ’ 

(2.30) 

with the projection operators 


(Ps)(Qia2)(/3i/32) ~ 2 A ^ai/32^a20i) ) 

(2.31) 

(Pa)(ai«2)(/ti/32) ~ 2 *^“2/32 ~ <^oi/32'^a2/3i) ) 

(2.32) 

which decompose the direct product space of two fundamental SU{Nc) 
tions into the irreducible representations 

represent a- 

N, ® Ne = (N, + l)Ne/2 © N,(N, - l)/2 . 

(2.33) 


With TrN^ 0 Nc = A^c the projector properties 

P^a ~ Ps,a) TrN^ 0 N^P 5 = (W+l)A^c/2, and TrN^ 0 NcPa = (A^c ~ l)A^c/2 , (2.34) 

we hnd for the loop-loop correlation function with both loops in the fundamental 
SU{Nc) representation 

(»'n.[Ci]»'n.[ 0])^ = exp + XS,&) (2.35) 



^Note that the Euclidean xSiSi 7 ^ 0 in contrast to xSiSi = 0 for Minkowskian light-like loops 
Ci considered in the original version of the Berger-Nachtmann approach [8,45]. 
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where 


(2.36) 


6*2 (N^ 


Nl-l 

2N, 


For one Wegner-Wilson loop in the fundamental and one in the adjoint represen¬ 
tation of SU{Nc), r'l = Nc and r 2 = N^ —1, which is needed in Sec. 4 to investigate 
the chromo-held distributions around color sources in the adjoint representation, the 
decomposition (2.29) reads 

TV 1 1 

® Cm = - Y + 2 “ 2 

with the projection operators^ Pi, P 2 , and P 3 that decompose the direct product 
space of one fundamental and one adjoint representation of SU{Nc) into the irre¬ 
ducible representations 

Ne®N2-l = Ne © ^N,(Ne-l)(Ne + 2) © 1 N,(Ne + 1)(N^ - 2), (2.38) 

which reduces for TV^ = 3 to the well-known SU (3) decomposition 


3©8 = 3©15©6. 


(2.39) 


With = Nc{Nf — 1) and projector properties analogous to (2.34), 

we obtain the loop-loop correlation function for one loop in the fundamental and 
one in the adjoint representation of SU{Nf) 


W^ACi]W^l-i[C 2 ]) =exp 
G 


CM 

n XS'iS’i © n \S2S2 


(2.40) 


X 


iV 2 -l 


exp 


TV, 


, iVc + 2 

+ 2(]VmI)“P 


2VS1S2 


+ 


iVe -2 

2(W-1) 


exp 


\_-jXS1S2 


where 

C'2(N2-1) = W. (2.41) 


Note that our expressions for the loop-loop correlation function (2.29) and, more 
specihcally, (2.35) and (2.40), are rather general results—as is our result for the VEV 
of one Wegner-Wilson loop (2.14)—obtained directly from the color neutrality of the 
QCD vacuum and the Gaussian approximation in the gluon held strengths. The loop 
geometries, which characterize the problem under investigation, are again encoded 

^The explicit form of the projection operators Pi, P 2 , and P 3 can be found in [51] but note that 
we use the Gell-Mann (conventional) normalization of the gluons. The eigenvalues, Aj, of the projec¬ 
tion operators in (2.37) can be evaluated conveniently with the computer program COLOUR [52]. 


11 



in the functions XSiSj, where also more detailed aspects of the QCD vacuum enter in 
terms of i.e., the gauge-invariant bilocal gluon held strength correlator (2.12). 

For the explicit computations of Xs’iS ’2 presented in Appendix C, one has to 
specify surfaces S' 1^2 with the restriction dSi ^2 = C'lg according to the non-Abelian 
Stokes theorem. We choose for S' 1^2 minimal surfaces that are built from the plane 
areas spanned by the corresponding loops Ci ^2 and the inhnitesimally thin tube 
which connects the two surfaces Si and 82 - This is in line with our surface choice in 
applications of the LLCM to high-energy reactions [8,32-34]. The thin tube allows 
us to compare the held strengths in surface Si with the held strengths in surface 82 - 

Due to the Gaussian approximation and the associated truncation of the cumu- 
lant expansion, the non-perturbative conhning contribution to the loop-loop cor¬ 
relation function depends on the surface choice. For example, our results for the 
chromo-held distributions of color dipoles obtained with the minimal surfaces dif¬ 
fer quantitatively from the ones obtained with the pyramid mantle choice for the 
surfaces [25] even if the same parameters are used. The qualitative main features 
of the non-perturbative SVM component (such as conhnement via hux tube forma¬ 
tion), however, emerge very similarly in both scenarios. From a comparison of the 
static quark-antiquark potential to the energy stored in the chromo-electric helds 
presented in Sec. 5, we infer that the minimal surfaces are more compatible with the 
Gaussian approximation. Indeed, the application of low-energy theorems in Sec. 5 
will show that the minimal surfaces are important for the consistency between the 
results for the VEV of one loop, (1F].[G]), and the loop-loop correlation function, 
{Wrf[Ci\Wr 2 [C 2 \) ■ In addition, the simplicity of the minimal surfaces gives dehnitive 
advantages in analytical computations. For example, it has allowed us to represent 
the conhning string as an integral over stringless dipoles with a given dipole number 
density [33]. 

In applications of the model to high-energy scattering [8,32-34] the surfaces are 
interpreted as the world-sheets of the conhning QGD strings in line with the picture 
obtained for the static dipole potential from the VEV of one loop. The minimal 
surfaces are the most natural choice to examine the scattering of two rigid strings 
without any huctuations or excitations. Our model does not choose the surface 
dynamically and, thus, cannot describe string hips between two non-perturbative 
color dipoles. Recently, new developments toward a dynamical surface choice and a 
theory for the dynamics of the conhning strings have been reported [53]. 
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2.3 Perturbative and Non-Perturbative QCD Components 


We decompose the gauge-invariant bilocal gluon field strength correlator (2.12)—as 
in the Minkowskian version of our model [8]—into a perturbative (P) and non- 
perturbative {NP) component 


p _ fpP pNP 

fii/per pupa ' ■‘■pupa ) 


(2.42) 


where gives the low-frequency background field contribution modeled by the 

non-perturbative stochastic vacuum model [21] and i'h® additional high-frequency 
contribution described by perturbative gluon exchange. This combination allows us 
to describe long and short distance correlations in agreement with lattice calcula¬ 
tions of the gluon field strength correlator [18,22-24]. Moreover, this two component 
ansatz leads to the static quark-antiquark potential with color Coulomb behavior 
for small and confining linear rise for large source separations in good agreement 
with lattice data as shown in Sec. 3. Note that in addition to our two component 
ansatz an ongoing effort to reconcile the non-perturbative SVM with perturbative 
gluon exchange has led to complementary methods [53-55]. 

We compute the perturbative correlator P^p^. from the Euclidean gluon propa¬ 
gator in the Feynman-’t Hooft gauge: 


(6;(Vi)ej(A'2)) 


-iK{Xi-X2) 

(27r)^ K'^ -I- uIq 


(2.43) 


where we introduce an effective gluon mass of mo = rOp = 0.77 GeV to limit the 
range of the perturbative interaction in the infrared (IR) region. This IR cutoff for 
the perturbative component is important in applications of our model to high-energy 
scattering [8]. Its value has been chosen such that the unintegrated gluon distribu¬ 
tion for transverse momenta below |A;_l| ~ 1 GeV is dominated by non-perturbative 
physics [33]. Of course, the parameter me is also important for the interplay be¬ 
tween the perturbative and non-perturbative components in the presented Euclidean 
applications. Furthermore, our value for mo gives the “perturbative glueball” {GB) 
generated by our perturbative component a finite mass of Mq^ = 2mc = 1.54 GeV, 
which is larger than that of its non-perturbative counterpart discussed below. This 
ensures that long-range correlations are dominated by non-perturbative physics. 

In leading order in the strong coupling g, the resulting bilocal gluon field strength 
correlator is gauge invariant already without the parallel transport to a common 


13 



reference point so that depends only on the difference Z = Xi — X 2 : 




(Z) = 


7r2 2 IdZ, 
'd^K 


{ZJ^, - Zp5^„) + ^ {Zp5,„ - Z^d.p)] Dp{Z^) (2.44) 


dZ^ ^ " 


-iKZ 


TT^J ip.TlY 


K,,KJ„ 


'■u-^'-cr'-' HP Ki,Kp5pa KpKpSycr KpK(^5yp 


b'p{K^ 


with the perturbative correlation function 

Dp{ZY = ^^,K2{mG\Z\) (2.45) 

/d*Ze-flp(Z^) = - 5 ^ . (2.46) 


The perturbative gluon held strength correlator has also been considered at 
next-to-leading order, where the dependence of the correlator on both the renor¬ 
malization scale and the renormalization scheme becomes explicit and an additional 
tensor structure arises together with a path dependence of the correlator [56]. How¬ 
ever, cancellations of contributions from this additional tensor structure have been 
shown [55]. We refer to Sec. 3.3 of Ref. [50] for a more detailed discussion of this 
issue. 


We describe the perturbative correlations in our phenomenological applications 
only with the leading tensor structure (2.44) and take into account radiative correc¬ 
tions by replacing the constant coupling with the running coupling 


= ^TiasiZY 


A 8 tt^ 

(33-2]Vain[(Z-2 + M2)/A|cD] 


(2.47) 


in the dual step of the computation of the y function, where the Euclidean dis¬ 
tance jZj over which the correlation occurs provides the renormalization scale. In 
Eq. (2.47) Nf denotes the number of dynamical quark havors, which is set to Nf = 0 
in agreement with the quenched approximation, Aqcd = 0.25 GeV, and M allows 
us to freeze g^ for jZj 00 . Relying on low-energy theorems, we freeze the running 
coupling at the value g"^ = 10.2 (= = 0.81), i.e. M = 0.488 GeV, at which our 

non-perturbative results for the confining potential and the total hux tube energy 
of a static quark-antiquark pair coincide (see Sec. 5). 

The tensor structure (2.44) together with the perturbative correlation func¬ 
tion (2.45) or (2.46) leads to the color Yukawa potential (which reduces for me = 0 
to the color Goulomb potential) as shown in Sec. 3. The perturbative contribution 
thus dominates the full potential at small quark-antiquark separations. 
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If the path connecting the points Xi and X 2 is a straight line, the non-perturbative 
correlator also depends only on the difference Z = Xi — X 2 . Then, the most 

general form of the correlator that respects translational, Lorentz, and parity invari¬ 
ance reads [21] 


pNP 


fiiypa 


(z) = fz:{z) + fzt{z) 


where 


(2.48) 


3(iV2 - 1 ) 


G 2 K D{Z^) 


+ (1 - ^ {ZJ^p - Z,6Z + {ZpS,^ - ZJZ] Di{Z^)] 

d^K 


d 


3(iV2-l)^^7 (27r) 
- (1 - k) 


^-iKZ 




K K 5 - K K 5 + K K 5 - K K 5 




D[{K^) := 


dK^ 


d'^ZDi{Z^)e 


2\ AKZ 


(2.49) 


In all previous applications of the SVM, this form, depending only on Z = Xi — X 2 , 
has been used. New lattice results on the path dependence of the correlator [57] 
show a dominance of the shortest path. This result is effectively incorporated in the 
model since the straight paths dominate in the averaging over all paths. 


Let us emphasize that the non-perturbative correlator (2.48) is a sum of the 
two different tensor structures, FZpa‘^ and FZpai with characteristic behavior: The 
tensor structure FZpa'^ is characteristic for Abelian gauge theories, exhibits the 
same tensor structure as the perturbative correlator (2.44) and does not lead to 
confinement [21]. In contrast, the tensor structure FZpa can occur only in non- 
Abelian gauge theories and Abelian gauge theories with monopoles and leads to 
confinement [21]. Therefore, we call the tensor structure multiplied by (1 — k) non¬ 
confining (nc) and the one multiplied by k confining (c). 

The non-perturbative correlator (2.48) involves the gluon condensate G 2 '■ = 

2 

[58], the weight parameter k, and the correlation length a which 
enters through the non-perturbative correlation functions D and Di. While the 
perturbative correlation function Dp given in (2.45) is computed from the gluon 
propagator (with a finite effective gluon mass), the non-perturbative correlation 
functions D and Di can be studied rigorously in lattice QCD investigations [18,22- 
24]. In addition, the non-perturbative correlation functions are constrained by the 
following physical considerations, (i) The correlations at large distances should de¬ 
crease exponentially so that the interaction range is determined by the glueball mass, 
(ii) Toward small distances, the non-perturbative correlation functions must satisfy 
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D(0) = -Di(O) = 1 to ensure the correct relation between the VEV of inhnitesimal 
plaquettes and the gluon condensate G 2 . (iii) The correlation functions must stay 
positive at all distances to be compatible with a spectral representation [59]. 

We adopt for our calculations a simple exponential correlation function 

D{Z^) = Di{Z^) = exp(-|Z|/a) , (2.50) 

which is consistent with the physical constraints discussed and has been successfully 
tested in hts to lattice data of the gluon held strength correlator [18,24]. The 
exponential correlation function stays positive for all Euclidean distances Z and 
is compatible with a spectral representation of the correlation function [59]. This 
is a conceptual improvement since the correlation function used in several earlier 
applications of the SVM becomes negative at large distances [12,25,31,40-46]. 

With the exponential correlation function (2.50), hts to the lattice data of the 
gluon held strength correlator down to distances of 0.4 fm give the following values 
for the parameters of the non-perturbative correlator [24]: G 2 = 0.173 GeV^, k = 
0.746, and a = 0.219 fm. We have optimized these parameters in a ht to high-energy 
scattering data® [ 8 ]: 


a = 0.302 fm, n = 0.74, G 2 = 0.074 GeV^ . (2.51) 

We use these optimized parameters (2.51) throughout this work. They lead to a 
static quark-antiquark potential that is in good agreement with lattice data (see 
Sec. 3) and, in particular, give a QGD string tension (3.12) of 0-3 = 0.22 GeV^ = 
1.12GeV/fm which is consistent with hadron spectroscopy [62], Regge theory [63], 
and lattice QGD investigations [64]. Moreover, the non-perturbative component 
with a = 0.302 fm generates a “non-perturbative glueball” with a mass of = 
2/a = 1.31 GeV which is smaller than Mq^ = 1.54 GeV and thus governs the long- 
range correlations as expected. We thus have one model that describes both static 
hadronic properties and high-energy reactions of hadrons and photons in good agree¬ 
ment with experimental and lattice QGD data. 

At this point, we would like to comment on the model parameters and the ac¬ 
curacy of the results. Although there are strong hints for the choice of the integra¬ 
tion surface and physical constraints on the non-perturbative correlation functions, 

^ Since we describe both lattice QCD data obtained in the quenched approximation and high- 
energy scattering data taken in the presence of light quarks, our value for the gluon condensate, 
G 2 = 0.074 GeV^, interpolates between G 2 = 0.173 GeV^ found in quenched lattice QCD inves¬ 
tigations [22,24] and G 2 = 0.024 ± 0.011 GeV^ found in phenomenology [58,60] and full lattice 
QCD investigations [23,24]. It is known that the effect of light quarks reduces the value of G 2 
substantially [61]. 
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Table 1; Different sets of parameters used together with different surfaces and different 
correlation functions in applications of the SVM to high-energy scattering. 


Reference 

|12| 

[41] 

[43] 

[8] 

a (fm) 

0.350 

0.346 

0.346 

0.302 

G 2 (GeV^) 

0.0605 

0.0631 

0.0631 

0.074 

K 

0.74 

0.74 

0.74 

0.74 

me (GeV) 

— 

— 

0.571 

0.77 

Integration surface 

pyramid 

pyramid 

pyramid 

minimal 

Gorrelation function 

Bessel 

Bessel 

Bessel 

exponential 

Perturbative component 

no 

no 

yes 

yes 


we have no criteria from first principles that £x these model ingredients unam- 
bigiously. Therefore, we have checked different integration surfaces and different 
non-perturbative correlation functions: While the analytic result for the string ten¬ 
sion changes, the general picture (e.g. the conhning linear rise of the static dipole 
potential and flux tube formation) is reproduced by readjusting the parameters a, 
K, and G 2 . Therefore, the model parameters are meaningful only within about 20% 
accuracy, which estimates possible errors incurred by chosing a certain combination 
of integration surfaces and correlation functions. In Table 1 we show different sets 
of parameters used together with different surfaces and different correlation func¬ 
tions in applications of the SVM to high-energy scattering. The table documents 
the stability of the SVM parameters within 20%. However, after the Gaussian ap¬ 
proximation (or truncation of the cumulant expansion) and the specihcation of the 
integration surface and the correlation functions, the quantitative results depend 
sensitively on some of the model parameters. To achieve a good £t to high-energy 
cross sections [8], a fine tuning of a and G 2 is necessary. 

Finally, we should discuss the pragmatic treatment of renormalization of the 
perturbative component (2.44) that dominates the small distance correlations. Only 
the lowest order result is adopted in which the strong coupling is promoted to a 1- 
loop effective running coupling. The mass renormalization of the considered heavy 
quarks and antiquarks is also taken into account by subtracting the self-energy 
of the sources in the computation of the static color dipole potential. Although 
phenomenologically successful, one needs to rehne the treatment of renormalization, 
for example, by explicitly taking into account counterterms for the cusps and by 
introducing a factorization scale in order to put the model on a more solid basis. 
We defer this task to future work and turn now to the phenomenological performance 
of our pragmatic approach. 
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3 The Static Color Dipole Potential 


In this section the QCD potential of static color dipoles in the fundamental and 
adjoint representation of SU{Nc) is computed in our model. Color Coulomb behavior 
is found for small dipole sizes and the conhning linear rise for large dipole sizes. 
Casimir scaling is obtained in agreement with lattice QCD investigations. 

The static color dipole—two static color sources separated by a distance R in 
a net color singlet state—is described by a Wegner-Wilson loop hbrlC] with a rect¬ 
angular path C of spatial extension R and temporal extension T ^ oo where r 
indicates the SU{Nc) representation of the sources considered. Figure 1 illustrates 
a static color dipole in the fundamental representation r = Nc. The potential of the 
static color dipole is obtained from the VEV of the corresponding Wegner-Wilson 
loop [15, 65] 

Vr{R) = - lim i ln(Wc[C'])pot , (3.1) 

T—>oo ± 

where “pot” indicates the subtraction of the self-energy of the color sources. This 
subtraction corresponds to the mass renormalization of the heavy color sources as 
discussed above. The static quark-antiquark potential is obtained from a loop 
in the fundamental representation (r = Nc) and the potential of a static gluino pair 
I/n2-i from a loop in the adjoint representation (r = N^ —1). 

With our result for (IW[C']), (2.14), obtained with the Gaussian approximation 



Figure 1: A static color dipole of size R in the fundamental representation. The rectan¬ 
gular path C of spatial extension R and temporal extension T indicates the world-line of 
the dipole described the Wegner-Wilson loop The shaded area bounded by the 

loop C = dS represents the minimal surface S used to compute the static dipole potential. 
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in the gluon field strength, the static potential reads 

yAR) = lim ^Xsspot (3.2) 

with the self-energy subtracted, i.e. Xsspot := Xss — Xssseif (see Appendix C). 
According to the structure of the gluon field strength correlator, (2.12) and (2.42), 
there are perturbative (P) and non-perturbative (AP) contributions to the static 
potential: 

K(P) = ^1^ ^ {xfspot + {XssAA + Xss At) ] ’ (3-3) 

where the explicit form of the y functions is given in (C.9), (C.28), and (C.37). 

The perturbative contribution to the static potential describes the color Yukawa 
potential (which reduces to the color Coulomb potential [66] for me = 0) 

Vf{R) = -C2(r) ^^expf-meP] . (3.4) 

Here we have used the result for Xss pot given in (C.37) and the perturbative corre¬ 
lation function 

:= f (3.5) 

J (27r)'’ 4t^\z\ 

which is obtained from the massive gluon propagator (2.43). As shown below, the 
perturbative contribution dominates the static potential for small dipole sizes R. 

The non-perturbative contributions to the static potential, the non-eonfining 
component (nc) and the confining component (c), read 

= C;(r) AtP-p D'Ym . (3.6) 

= g;(r) j\(R-p)D^^\p-‘) . (3.7) 

where we have used the results for x^spil Xsfpot = Xss given respectively 
in (C.28) and (C.9) obtained with the minimal surface, i.e. the planar surface 
bounded by the loop as indicated by the shaded area in Fig. 1. With the expo¬ 
nential correlation function (2.50), the correlation functions in (3.6) and (3.7) read 

D'P\Z^) := j = - o |ZpA',||Z|/o] , (3.8) 

D»>(Z=) := j D(K^)MKi) = 2 |Z| A'i[|i|/al . (3.9) 
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For large dipole sizes, R > 0.5 fm, the non-confining contribution (3.6) vanishes 
exponentially while the confining contribution (3.7)—as anticipated—leads to con¬ 
finement [21], i.e. the confining linear increase. 


\R) 


ij>0.5fm 


= arR + const. 


(3.10) 


Thus, the QCD string tension is given by the confining SVM component [21]: For 
a color dipole in the SU{Nf) representation r, it reads 


(T Y" 


C2 {t) 


TT^G2K 

48 



dZ^D{Z^) 


C2(r) 


7r^KG2a‘^ 

24 


(3.11) 


where the exponential correlation function (2.50) is used in the final step. Since 
the string tension can be computed from first principles within lattice QCD [64], 
relation (3.11) puts an important constraint on the three parameters of the non- 
perturbative QCD vacuum a, G 2 , and k. With the values for a, G 2 , and k given 
in (2.51), which are used throughout this work, one obtains for the string tension of 
the SU (3) quark-antiquark potential (r = 3) a reasonable value of 


as = 0.22 GeV^ = 1.12 GeV/fm . 


(3.12) 


The static SU{Nc = 3) quark-antiquark potential Vnc(-R) = VsiR) is shown as 
a function of the quark-antiquark separation R in Fig. 2, where the solid, dotted, 
and dashed lines indicate the full static potential and its perturbative and non- 
perturbative contributions, respectively. For small quark-antiquark separations R < 
0.5 fm, the perturbative contribution dominates giving rise to the well-known color 
Coulomb behavior. For medium and large quark-antiquark separations R > 0.5 fm, 
the non-perturbative contribution dominates and leads to the confining linear rise 
of the static potential. The transition from perturbative to string behavior takes 
place at source separations of about 0.5 fm in agreement with the recent results 
of Liischer and Weisz [26]. This supports our value for the gluon mass me = 
mp = 0.77 GeV which is important only around R ~ 0.4 fm, i.e. for the interplay 
between perturbative and non-perturbative physics. For R < 0.3 fm and R > 0.5 fm, 
the effect of the gluon mass, introduced as an IR regulator in our perturbative 
component, is negligible. String breaking is expected to stop the linear increase 
for i? > 1 fm where lattice investigations show deviations from the linear rise in 
full QCD [64,67]. As our model is working in the quenched approximation, string 
breaking through dynamical quark-antiquark production is excluded. 

As can be seen from (3.2), the static potential shows Casimir scaling which 
emerges in our approach as a trivial consequence of the Gaussian approximation 
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0 0.5 1 1.5 


R [fm] 

Figure 2: The static SU{Nc = 3) quark-antiquark potential V^^{R) = V 3 {R) as a function 
of the quark-antiquark separation R. The solid, dotted, and dashed lines indicate the full 
static potential and its perturbative and non-perturbative contributions, respectively. For 
small quark-antiquark separations, R < 0.5 fm, the perturbative contribution dominates 
and gives rise to the well-known color Coulomb behavior at small distances. For medium 
and large quark-antiquark separations, R > 0.5 fm, the non-perturbative contribution 
dominates and leads to the confining linear rise of the static potential. As our model is 
working in the quenched approximation, string breaking cannot be described, which is 
expected to stop the linear increase for i? > Ifm [64,67]. 


used to truncate the cumulant expansion (2.7). Indeed, the Casimir scaling hypoth¬ 
esis [68] has been verihed to high accuracy for SU{3) on the lattice [19,20] (see also 
Fig. 3). These lattice results have been interpreted as a strong hint toward Gaus¬ 
sian dominance in the QCD vacuum and thus as evidence for a strong suppression 
of higher cumulant contributions [69,70]. In contrast to our model, the instanton 
model can describe neither Casimir scaling [70] nor the linear rise of the confining 
potential [71]. 

Figure 3 shows the static SU{Nc = 3) potential for fundamental sources VN^(i?) = 
V^{R) (solid line) and adjoint sources lN 2 _i(i?) = Vs{R) (dashed line) as a function 
of the dipole size R in comparison to SU{3) lattice data [20,64]. The model results 
are in good agreement with the lattice data. In particular, the obtained Casimir 
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Figure 3: The static SU{Nc = 3) potential of color dipoles in the fundamental represen¬ 
tation V 3 (i?) (solid line) and adjoint representation Vs{R) (dashed line) as a function of 
the dipole size R in comparison to SU{3) lattice data for (3 = 6.0, 6.2, and 6.4 [20,64]. 
The model results are in good agreement with the lattice data. This particularly holds for 
the obtained Casimir scaling behavior. 


scaling behavior is strongly supported by SU{3) lattice data [19,20]. This, however, 
points also to a shortcoming of our model: From Eq. (3.2) and Fig. 3 it is clear that 
string breaking is described neither for fundamental nor for adjoint dipoles in our 
model which indicates that not only dynamical fermions (quenched approximation) 
but also some gluon dynamics are missing. 

An extension of the model that allows one to describe color screening remains a 
major challenge. Without such a modihcation, our model unfortunately cannot con¬ 
tribute to the recent discussion on the scaling behavior of fc-string tensions i.e. 
the tensions of strings connecting sources with A/'-ality A; > 1. A source of A/'-ality 
k < Nc/2 is dehned as a source in the representation constructed from the tensor 
product of quarks—objects transforming under the fundamental representation— 
and antiquarks—objects transforming under the conjugated representation—where 
k is the number of quarks minus the number of antiquarks modulo N^, see e.g. [72]. 
For SU{Nc) with Nc > 4:, k strings are particularly interesting since in addition 
to the fundamental {k = 1) string other strings also exist that are stable against 
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color screening. Based on lattice results for Nc = 4,5,6 [72,73], the present de¬ 
bate is whether the corresponding string tensions show Casimir scaling behavior 
(Tfc oc k{Nc — k) , the sine law behavior ak oc sm.{'Kk/Nc) predicted from M-theory 
approaches to QCD [74], or simply the behavior of k non-interacting fundamental 
strings ak = ka-^^. Physical explanations of the lattice results obtained are discussed, 
for example, in the center vortex conhnement mechanism [7, 75]. 


4 Chromo-Field Distributions of Color Dipoles 


In this section we compute the chromo-electric helds generated by a static color 
dipole in the fundamental and adjoint representation of SU{Nc). We hnd formation 
of a color flux tube that conhnes the two color sources in the dipole. This conhning 
string is analyzed quantitatively. Its mean squared radius is calculated and trans¬ 
verse and longitudinal energy density prohles are provided. The interplay between 
perturbative and non-perturbative contributions to the chromo-field distributions is 
investigated and exact Casimir scaling is found for both contributions. 

As already explained in Sec. 3, the static color dipole—two static color sources 
separated by a distance 7? in a net color singlet state—is described by a Wegner- 
Wilson loop ITr[C'] with a rectangular path C of spatial extension R and temporal 
extension T —cx) (cf. Fig. 1) where r indicates the SU{Nc) representation of the 
sources considered. A second small quadratic loop or plaquette in the fundamental 
representation placed at the space-time point X with side length Rp —0 and 
oriented along the af3 axes. 


= TrN.exp 


-w 


dZ^G' 


\{Z)t 


Nc 


Cp 


= l-R\ 


AN, 


+ 0{R 


(4.1) 


is needed—as a “Hall probe”—to calculate the chromo-held distributions at the 
space-time point X caused by the static sources [76, 77] 
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{UV[C]P™(A)) 
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N"f(A')> 


(4.2) 

(4.3) 


with no summation over a and (3 in (4.1), (4.2), and (4.3). In dehnition (4.2) 
(• • ■)wr[c] indicates the VEV in the presence of the static color dipole while (.. .)vac 
indicates the VEV in the absence of any color sources. Depending on the plaquette 


23 



orientation indicated by a and (3, one obtains from (4.3) the squared components of 
the chromo-electric and chromo-magnetic helds at the space-time point X: 
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(4.4) 


i.e. space-time plaquettes (a/? = i4) measure chromo-electric helds and space-space 
plaquettes (a/? = ij) chromo-magnetic helds. As shown in Fig. 4, we place the static 
color sources on the Xi axis at (Xi = ±i?/2, 0, 0, X4) and use the following notation 
plausible from symmetry arguments 


El = El = El, Bi = = Bj . (4.5) 

Figure 4 illustrates also the plaquette at X = (Xi,X2,0,0) needed to com¬ 

pute P|^(X). Due to symmetry arguments, the complete information on the chromo- 
held distributions is obtained from plaquettes in “transverse” space X = (Xi, X2, 0, 0) 
with four diherent orientations, af3 = 14, 24, 13, 23 [cf. Eq. (4.5)]. 

The energy and action density distributions around a static color dipole in the 
SU (Xc) representation r are given by the squared chromo-held distributions 


SriX) 

Sr{X) 


-E^{X) + B\X) 
■ \e‘^{X) + B^{X) 


(4.6) 

(4.7) 


with signs according to Euclidean space-time conventions. Low-energy theorems 
that relate the energy and action stored in the chromo-helds of the static color 
dipole to the corresponding ground state energy are discussed in the next section. 

For the chromo-held distributions of a static color dipole in the fundamental 
representation of SU{Nc), i.e. a static quark-antiquark pair, we obtain with our 
results for the VEV of one loop (2.14) and the correlation of two loops in the 
fundamental representation (2.35) 
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Figure 4: The plaquette-loop geometry needed to compute the squared chromo-electric 
field -Kji (X) generated by a static color dipole in the fundamental SU{Nc) representation 
(r = Nc). The rectangular path C indicates the world-line of the static dipole described 
the Wegner-Wilson loop VFnc[C']. The square with side length Rp illustrates the plaquette 
Pj^^(X). The shaded areas represent the minimal surfaces used in our computation of the 
chromo-field distributions. The thin tube allows us to compare the gluon field strengths 
in surface Sp with the gluon field strengths in surface Sw- 

where XSiSj is defined in (2.27). The subscripts P and W indicate surface integra¬ 
tions to be performed over the surfaces spanned by the plaquette and the Wegner- 
Wilson loop, respectively. Choosing the surfaces—as illustrated by the shaded areas 
in Fig. 4—to be the minimal surfaces connected by an infinitesimal thin tube (which 
gives no contribution to the integrals) it is clear that XSpSp R% and XSpSw ^p- 
Being interested in the chromo-fields at the space-time point X, the extension of 
the quadratic plaquette is taken to be infinitesimally small, Rp —> 0, so that one 
can expand the exponential functions and keep only the term of lowest order in Rp: 

= -C,m T ^ . (4.9) 

This result—obtained with the matrix cumulant expansion in a very straightforward 
way—agrees exactly with the result derived in [25] with the expansion method. 
Indeed, the expansion method agrees for small y functions with the matrix cumulant 
expansion (Berger-Nachtmann approach) used in this work but breaks down for large 
y functions, where the matrix cumulant expansion is still applicable. 


25 
















The chromo-field distributions of a static color dipole in the adjoint represen¬ 
tation of SU{Nc), i.e. a static gluino pair, are computed analogously. Using our 
result (2.40) for the correlation of one loop in the fundamental representation (pla- 
quette) with one loop in the adjoint representation (static sources), one obtains 
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which reduces—as explained for sources in the fundamental representation—to 




(4.11) 


Thus, the squared chromo-electric fields of an adjoint dipole differ from those of a 
fundamental dipole only in the eigenvalue of the corresponding quadratic Casimir 
operator C 2 (rj. In fact, Casimir scaling of the chromo-field distributions holds for 
dipoles in any representation r of SU (iVc) in our model. As can be seen with the low- 
energy theorems discussed below, this is in line with the Casimir scaling of the static 
dipole potential found in the previous section. In addition to lattice investigations 
that show Casimir scaling of the static dipole potential to high accuracy in SU (3) [19, 
20], Casimir scaling of the chromo-field distributions has been considered on the 
lattice as well but only for SU{2) [78]. Here only slight deviations from the Casimir 
scaling hypothesis have been found which were interpreted as hints toward adjoint 
quark screening. 


In our model the shape of the field distributions around the color dipole is iden¬ 
tical for all SU{Nc) representations r and given by XspSw' This again illustrates 
the shortcoming of our model discussed in the previous section. Working in the 
quenched approximation, one expects a difference between fundamental and adjoint 
dipoles: string breaking cannot occur in fundamental dipoles as dynamical quark- 
antiquark production is excluded but should be present for adjoint dipoles because 
of gluonic vacuum polarization. Comparing (4.9) with (4.11) it is clear that this dif¬ 
ference is not described in our model. In fact, as shown in Sec. 3, string breaking is 
described neither for fundamental nor for adjoint dipoles. Interestingly, even on the 
lattice there has been no striking evidence for adjoint quark screening in quenched 
QCD [79]. It is even conjectured that the Wegner-Wilson loop operator is not suited 
to studies of string breaking [80]. 

In the LLCM there are perturbative (P) and non-perturbative (AP) contri¬ 
butions to the chromo-electric fields according to the structure of the gluon field 
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(4.12) 


strength correlator, (2.12) and (2.42), 

X ([xf,s„(A)]h+{ [xTji(x)]^^ + . 

Interference of perturbative and non-perturbative correlations is not considered to 
be in line with the applications of our model to high-energy scattering [8, 32-34] with 
separate hard (perturbative) and soft (non-perturbative) Pomeron exchanges. The 
interferences do not change the qualitative picture. Slight modihcations occur in 
regions where helds originating from perturbative and non-perturbative correlations 
are of similar size. For the y functions in (4.12), we give directly in the following 
the hnal results obtained with the minimal surfaces shown in Fig. 4. Details of their 
derivation can be found in Appendix C. 

The perturbative contribution (P) described by massive gluon exchange leads, of 
course, to the well-known color Yukawa fieldwhich. reduces to the color Coulomb field 
for niG = 0. It contributes only to the chromo-electric helds, {afi = 14) 

and P^ = Ey = El {ajd = 24), and reads explicitly for X = (Xi, X 2 , 0, 0) 

d 2 poo 

= -^jiT{(Xi-Rmf(ZlA)Dp(ZlV 

-(X, + RI2)g\Zlc)Dp(Zlc)} ( 4 . 13 ) 

d 2 poo 

[XspsPX%_^ ^ -fij dTX.ig^Z'ii) Dp(Z'i^) - g'\Z'fc) DpiZ'ic)] ( 4 . 14 ) 
with the perturbative correlation function (2.45), the running coupling (2.47), and 

ZIa = ) + ^2 + and = (^i + f) + ^2 + . (4.15) 


The non-confining non-perturbative contribution {NPnc) has the same structure 
as the perturbative contribution—as expected from the identical tensor structure— 
but differs, of course, in the prefactors and the correlation function, Di 7 ^ Dp. Its 
contributions to the chromo-electric helds Pj^ = El {afi = 14) and El_ = E^ = El 
(ck/ 5 = 24) read for X = (Xi, X 2 , 0, 0) 




flp>r“G2(l-K) 
6(iV|-l) , 


dT[(Xi - R/2) Di(Zlt) 


- (X, + fi/2)D,(Zfc)} (4.16) 


[xS(A')]24 = 


Rl7r^G^(l-n) 


dTXVD,{ZiA-D,(Z{c)} (4.17) 
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with the exponential correlation function (2.50) and Z\j^ and as given in (4.15). 

The confining non-perturbative contribution {NP c) has a different structure that 
leads to conhnement and flux-tube formation. It gives contributions only to the 
chromo-electric held {aj3 = 14) which read for X = (Xi, X 2 , 0, 0) 

[xS.(A')].4 = lypD'-^HZi) , (4.18) 

with the correlation function given in (3.9) as derived from the exponential correla¬ 
tion function (2.50), and 


Zi = [Xi + (1/2 - p)Rf + X| . (4.19) 


In our model there are no contributions to the chromo-magnetic fields, i.e. the 
static color charges do not affect the magnetic background held 

Bf^= Bl = 0 and Bl = Bl = = 0 , (4.20) 

which can be seen from the corresponding plaquette-loop geometries as pointed out 
in Appendix C. Thus, the energy and action densities are identical in our approach 
and completely determined by the squared chromo-electric helds 

efiX) = Sr{X) = ~E^{X) . (4.21) 

This picture is in agreement with other ehective theories of conhnement such as the 
’t Hooft-Mandelstam picture [81] or dual QCD [82] and, indeed, a relation between 
the dual Abelian Higgs model and the SVM has been established [83]. In con¬ 
trast, lattice investigations work at scales at which the chromo-electric and chromo- 
magnetic helds are of similar magnitude [30,84-86]. Indeed, these simulations have 
been performed in SU{Nc = 2) at bare couplings Qq down to fio = A/= 2.74 corre¬ 
sponding to a minimum lattice cutoh of 0.04 fm, which determines also the minimum 
size of the plaquette used in the measurements of the color helds. Interestingly, as 
shown in Figs. 13 and 14 of [84], the action density slightly decreases by decreasing 
the lattice spacing from 0.08 fm {Pq = 2.5) to 0.05 fm {Pq = 2.635). In our model the 
vanishing of the chromo-magnetic helds determines the value of the Callan-Symanzik 
P function at the renormalization scale at which our non-perturbative component is 
working. This is shown as a result of low-energy theorems in the next section. 

In Fig. 5 the energy density distributions g'^efiXi, X 2 = Xfi) generated by a color 
dipole in the fundamental SU{Nc = 3) representation (r = 3) are shown for quark- 
antiquark separations of R = 0.1, 0.5, 1 and 1.5 fm. With increasing dipole size 
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R = 0.1 fm 


R = 1 fm 



Figure 5: Energy density distributions g^ez{Xi,X 2 = X 3 ) generated by a color dipole 
in the fundamental SU{3) representation (r = 3) for quark-antiquark separations of 
R = 0.1, 0.5, 1 and 1.5 fm. Flux-tube formation leads to the confining QCD string with 
increasing dipole size R. 

R, one sees explicitly the formation of the flux tube which represents the conhning 
QCD string. 

The longitudinal and transverse energy density profiles generated by a color 
dipole in the fundamental representation (r = 3) of SU{Nc = 3) are shown for 
quark-antiquark separations (dipole sizes) of R = 0.1, 0.5, 1 and 1.5 fm in Figs. 6 
and 7. The perturbative and non-perturbative contributions are given by the dotted 
and dashed lines, respectively, and the sum of both in the solid lines. The open and 
hlled circles indicate the quark and antiquark positions. As can be seen from (4.3) 
and (4.4), we cannot compute the energy density separately but only the product 
g‘^er{X). Nevertheless, a comparison of the total energy stored in chromo-electric 
helds to the ground state energy of the color dipole via low-energy theorems yields 
g‘^ = 10.2 (= as = 0.81) for the non-perturbative SVM component as shown in the 
next section. 

In Figs. 6 and 7 the formation of the conhning string (hux tube) with increasing 
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Figure 6: Longitudinal energy density profiles g‘^e 3 {Xi,X 2 = X 3 = 0) generated by a color 
dipole in the fundamental SU (3) representation (r = 3) for quark-antiquark separations 
of i? = 0.1, 0.5, 1 and 1.5 fm. The dotted and dashed lines give the perturbative and 
non-perturbative contributions, respectively, and the solid lines the sum of both. The 
open and filled circles indicate the quark and antiquark positions. For small dipoles, R = 
0.1 fm, perturbative physics dominates and non-perturbative correlations are negligible. 
For large dipoles, R > 1 fm, the formation the confining string (flux tube) can be seen 
which dominates the chromo-electric fields between the color sources. 
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Figure 7: Transverse energy density profiles g‘^£ 3 {X 2 , Xi =^3 = 0 ) generated by a color 
dipole in the fundamental SU{3) representation (r = 3) for quark-antiquark separations of 
R = 0.1, 0.5, 1 and 1.5 fm. The dotted and dashed lines give the perturbative and non- 
perturbative contributions, respectively, and the solid lines the sum of both. The filled 
circles indicate the positions of the color sources. For small dipoles, R = 0.1 fm, perturba¬ 
tive physics dominates and non-perturbative correlations are negligible. For large dipoles, 
7? > 1 fm, the formation the confining string (flux tube) can be seen which dominates the 
chromo-electric fields between the color sources. 
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Figure 8: Root mean squared radius Rms of the flux tube and energy density in the 
center of a fundamental SU{3) dipole g^e^{X = 0) as a function of the dipole size R. 
Perturbative and non-perturbative contributions are given respectively by the dotted and 
dashed lines and the sum of both in the solid lines. For large R, both the width and 
height of the flux tube in the central region are governed completely by non-perturbative 
physics and saturate respectively at ~ 0.55 fm and = 0) IGeV/fm^. 

The latter value is extracted with the result = 10.2 deduced from low-energy theorems 
in the next section. 


source separations R can again be seen explicitly: For small dipoles, R = 0.1 fm, 
perturbative physics dominates and non-perturbative correlations are negligible. For 
large dipoles, i? > 1 fm, the non-perturbative correlations lead to formation of 
a narrow flux tube which dominates the chromo-electric helds between the color 
sources. 


Figure 8 shows the evolution of the transverse width (upper plot) and height 
(lower plot) of the flux tube in the central region of the Wegner-Wilson loop as a 
function of the dipole size R where perturbative and non-perturbative contributions 
are given by the dotted and dashed lines, respectively, and the sum of both in the 
solid lines. The width of the flux tube is best described by the root mean squared 
(ms) radius 


R 


ms 


f dX^Xlg^er(X, = 0,X^) 
f dX^X^g^er(X, = 0,X^) ’ 


(4.22) 
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which is universal for dipoles in all SU (iVc) representations r as the Casimir factors 
divide out. The height of the flux tube is given by the energy density in the center 
of the considered dipole, g‘^er{X = 0). For large source separations, i? > 1 fm, both 
the width and height of the flux tube in the central region of the Wegner-Wilson loop 
are governed completely by non-perturbative physics and saturate for a fundamental 
SU (3) dipole (r = Nc = 3) at reasonable values of 

~ 0.55 fm and = 0) ^ 1 GeV/fm^ with = 10.2 . (4.23) 

Note that the qualitative features of the non-perturbative SVM component do 
not depend on the specihc choice for the parameters, surfaces, and correlation func¬ 
tions and have already been discussed with the pyramid mantle choice of the surface 
and different correlation functions in the hrst investigation of flux-tube formation 
in the SVM [25]. The quantitative results, however, are sensitive to the parameter 
values, the surface choice, and the correlation functions and are presented above 
with the LLCM parameters, the minimal surfaces, and the exponential correlation 
function. 


5 Low-Energy Theorems 

In this section we use low-energy theorems to test the consistency of the non- 
perturbative SVM component and to determine the value of the Callan-Symanzik 
(3 function and = 9^! (47r) at the renormalization scale at which this component 
is working. The energy and action sum rules considered allow us to conhrm the 
consistency of our loop-loop correlation result with the result obtained for the VEV 
of one loop. Finally, we compare our results for (3 and as to model independent 
QCD results for the Callan-Symanzik (3 function. 

Many low-energy theorems have been derived in continuum theory by Novikov, 
Shifman, Vainshtein, and Zakharov [61] and in lattice gauge theory by Michael [27]. 
Here we consider the energy and action sum rules—known in lattice QCD as Michael 
sum rules —that relate the energy and action stored in the chromo-helds of a static 
color dipole to the corresponding ground state energy [15,65] 

= - hm lln{»V|C]> , (5.1) 

1 -^oo ± 

In their original form [27], however, the Michael sum rules are incomplete [28,31]. In 
particular, signihcant contributions to the energy sum rule from the trace anomaly 
of the energy-momentum tensor have been found [28] that modify the naively ex¬ 
pected relation in line with the importance of the trace anomaly found for hadron 
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masses [87]. Taking all these contributions into account, the energy and action sum 
rule read respectively [28-30] 


Er{R) = j 

Er(R) + R 


d^XeriX) 

dEr{R) _ 
OR ~ 



(5.2) 

(5.3) 


where the Callan-Symanzik function is denoted by j3{g) = ladg/dja with the renor¬ 
malization scale fi. 


Inserting (5.3) into (5.2), we hnd the following relation between the total energy 
stored in the chromo-helds i7);°*(i?) and the ground state energy Er{R): 

Er(R) := j SXEr(X) = i . (5.4) 

The difference from the naive classical expectation that the full ground state energy 
of the static color sources is stored in the chromo-helds is due to the trace anomaly 
contribution [28] described by the second term on the RHS of (5.2). Indeed, for 
the Coulomb potential, obtained in tree-level perturbation theory, the action sum 
rule (5.3) shows explicitly that the trace anomaly contribution vanishes on the clas¬ 
sical level as expected. 


With the low energy theorems (5.3) and (5.4) the ratio of the integrated squared 
chromo-magnetic to the integrated squared chromo-electric held distributions can 
be derived 

jd^XBHX) [2 + 3l3(g)lg]EAR) + [-i-P(g)l9]R^-Sgl 
fiPXmX) [2-3l3{g)/g]E4S) + l2 + l3(g)/g]R^' 

This ratio can be used, for example, to determine non-perturbatively the Callan- 
Symanzik j3{g) function. For SU{Nc = 2) lattice investigations along these lines 
have already been performed [30, 86,88, 89].® 

In the large R region, the static color dipole potential can be approximated by 
the linear potential Vr{R) = UrR = Er{R) — Egeu with string tension in the 
considered representation r. In this approximation, the ratio (5.5) becomes the 
simple form 

_ 2 -I- (3{g)/g 
Vr{R)=arR 2 - I3{g)/g ' 

®In [84] the /? function was determined similarly based on a high-statistics study of chromo-field 
distributions in SU{Nc = 2) but unfortunately without taking the trace anomaly contribution into 
account. 
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Since the non-perturbative SVM component of our model describes the conhning 
linear potential for large source separations R, we can use (5.6) together with the 
vanishing of the chromo-magnetic helds (4.20) to determine the value of the Callan- 
Symanzik (3 function at the scale at which the non-perturbative component is 
working 


Pi.g) 


9 




-2 . 


(5.7) 


Here one should emphasize that this value is strictly valid only at asymptotically 
large values of R, while perturbative correlations must be taken into account to 
extend this investigation to smaller values of R. 


Concentrating on the conhning non-perturbative component {NPc) we now use 
(5.4) to determine the value of a.^ = j (47r) at which the non-perturbative SVM 
component is working. The RHS of (5.4) is obtained directly from the conhning 
contribution to the static potential E^^^{R) = given in (3.7) in Sec. 3. The 

left-hand side (LHS) of (5.4), however, involves a division by the a priori unknown 
value of after integrating g‘^er{X) for the chromo-electric held of the conhning 
non-perturbative component (4.18). As discussed in the previous section, we cannot 
compute the energy density separately but only the product g‘^er{X). Adjusting the 
value of g^ such that (5.4) is exactly satished for source separations of R = 1.5 fm, 
we hud that the non-perturbative component is working at the scale /iArp at which 


g‘^{g-Np) — 10.2 = as{g-Np) — o.8i. 


(5.8) 


As already mentioned in Sec. 2.3, we use this value as a practical asymptotic limit for 
the simple one-loop coupling (2.47) used in our perturbative component. Note that 
earlier SVM investigations along these lines found a smaller value of as{fiNp) = 0.57 
with the pyramid mantle choice for the surface [25, 31] but were incomplete since 
only the contribution from the traceless part of the energy-momentum tensor was 
considered in the energy sum rule. 

In Fig. 9 we show the total energy stored in the chromo-held distributions around 
a static color dipole in the fundamental (r = 3) and adjoint (r = 8) representation 
of SU{3) from the conhning non-perturbative SVM component, for 

CKs = 0.81 (solid lines) as a function of the dipole size R. Comparing this total 
energy, which appears on the LHS of (5.4), with the corresponding RHS of (5.4) 
(dashed lines), we hud good consistency even down to very small values of R. This 
is a nontrivial and important result as it conhrms the consistency of our loop-loop 
correlation result—needed to compute the chromo-electric held—with the result 
obtained for the VEV of one loop—needed to compute the static potential 
Moreover, it shows that the minimal surfaces ensure the consistency of our non- 
perturbative component. The good consistency found for the pyramid mantle choice 
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Figure 9: The total energy stored in the chromo-field distributions around a static color 
dipole of size R in the fundamental (r = 3) and adjoint (r = 8) representation of SU (3) 
from the confining non-perturbative SVM component, for ag = 0.81 (solid 

lines) compared with the relation to the corresponding ground state energy (dashed lines) 
given by the low-energy theorem (5.4). Good consistency is found even down to very small 
values of R. 

of the surface relies on the naively expected energy sum rule [25, 31] in which the 
contribution from the traceless part of the energy-momentum tensor is not taken 
into account. 

Let us discuss the values oi P/g and at the renormalization scale ptvp— given 
respectively in (5.7) and (5.8)—in comparison with the perturbative expansion [90] 
and lattice computations [91] of the Callan-Symanzik function in pure SU{Nc = 3) 
gauge theory. We obtained (5.7) and (5.8) such that the renormalization scales 
fiNP appearing in the two equations should be in good agreement. Considering 
P/g as a. function g"^, one thus can compare our combination with the perturbative 
expansion [90]. This comparison shows that our result is close to the curve obtained 
on the two-loop level in perturbation theory. In contrast, the non-perturbative 
lattice results for the P function of Liischer et al. [91] are in good agreement with 
the perturbative three-loop result computed in the minimal modihed subtraction 
(MS) scheme [92], However, it must be stressed that in the lattice investigation the 
considered values of the running coupling g'^ stay below 3.5 while our comparison 
requires values up to g^{fi]\fp) = 10.2. Thus, relying on a large extrapolation of the 
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model independent QCD results, our comparison provides at best an orientation. For 
a meaningful consistency check, we have to map out the Callan-Symanzik function at 
smaller values of i?, where also perturbative correlations must be taken into account 
and thus rehnements of our treatment of renormalization are needed. The low-energy 
theorems will provide crucial criteria for the success of such improvements. 


6 Euclidean Approach to High-Energy Scattering 


In this section we present a Euclidean approach to high-energy reactions of color 
dipoles in the eikonal approximation. After a short review of the functional integral 
approach to high-energy dipole-dipole scattering in Minkowski space-time, we gen¬ 
eralize the analytic continuation introduced by Meggiolaro [35] from parton-parton 
scattering to dipole-dipole scattering. This shows how one can access high-energy 
reactions directly in lattice QCD. We apply this approach to compute the scattering 
of dipoles in the fundamental and adjoint representation of SU (W) at high-energy 
in the Euclidean LLCM. The result shows the consistency with the analytic contin¬ 
uation of the gluon held strength correlator used in all earlier applications of the 
SVM and LLCM to high-energy scattering. Finally, we comment on the QCD van 
der Waals potential which appears in the limiting case of two static color dipoles. 

In Minkowski space-time, high-energy reactions of color dipoles in the eikonal ap¬ 
proximation have been considered—as basis for hadron-hadron, photon-hadron, and 
photon-photon reactions—in the functional integral approach to high-energy colli¬ 
sions developed originally for parton-parton scattering [9,10] and then extended to 
gauge-invariant dipole-dipole scattering [11-13]. The corresponding T-matrix ele¬ 
ment for the elastic scattering of two color dipoles at transverse momentum transfer 
Q± (t = —q^) and c.m. energy squared s reads 


T^r 2 i^, t, Z 2 , f 2 i) = ‘2is 




l-S^ {s,h^,zi,ri^,Z2,f2±) ( 6 . 1 ) 


with the S'-matrix element (M refers to Minkowski space-time) 




M 


(s, h^,zi, ri_L, 2 : 2 , rax) = hm 

T—>00 


{Wr,[Cl]WrAC2])M 

{Wr,[Cl])M{Wr,[C2])M 


( 6 . 2 ) 


The color dipoles are considered in the SU (W) representation r* and have transverse 
size and orientation The longitudinal momentum fraction carried by the quark 
of dipole i is Zi. [Here and in the following we use several times the term quark 
generically for color sources in an arbitary SU{N^ representation.] The impact 
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(6.3) 


parameter between the dipoles is [44] 

h = fig + (1 - Zi)fi± - r 2 q - (1 - Z 2 )f 2 ± = ricm- ^ 2 cm , 

where r^g (rjg) is the transverse position of the quark (antiquark), = fig — 
fig, and ficm = Zifig + (1 — Zi)fiq is the center of light-cone momenta. Figure 10 
illustrates the (a) space-time and (b) transverse arrangement of the dipoles. The 
dipole trajectories Q are described as straight lines. This is a good approximation 
as long as the kinematical assumption behind the eikonal approximation, s 3> —t, 
holds which allows us to neglect the change of the dipole velocities Vi = Pi/m m. 
the scattering process, where Pi is the momentum and m the mass of the considered 
dipole. Moreover, the paths Q are considered light-like^ in line with the high-energy 
limit, w? s —> cx). For the hyperbolic angle or rapidity gap between the dipole 
trajectories 7 = (ui ■ V 2 )—which is the central quantity in the analytic continuation 
discussed below and also dehned through s = Anif cosh^( 7 / 2 )—the high-energy limit 
implies 

lim 7 ln(s/m^) —> 00 . (6.4) 

The QCD VEV’s (.. .)m in the F'-matrix element (6.2) represent Mm/cows/czan func¬ 
tional integrals [10] in which—as in the Euclidean case discussed above—the func¬ 
tional integration over the fermion helds has already been carried out. 

The Euclidean approach to the described elastic scattering of dipoles in the 
eikonal approximation is based on Meggiolaro’s analytic continuation of the high- 
energy parton-parton scattering amplitude [35]. Meggiolaro’s analytic continuation 
has been derived in the functional integral approach to high-energy collisions [9,10] 
in which parton-parton scattering is described in terms of Wegner-Wilson lines: 
The Minkowskian amplitude, g^{'-),T,t), given by the expectation value of two 
Wegner-Wilson lines, forming a hyperbolic angle 7 in Minkowski space-time, and the 
Euclidean “amplitude,” g^{Q,T,t), given by the expectation value of two Wegner- 
Wilson lines, forming an angle 0 G [0, tt] in Euclidean space-time, are connected 
by the following analytic continuation in the angular variables and the temporal 
extension T, which is needed as an IR regulator in the case of Wegner-Wilson lines, 

g^{Q,T,t) = g^{^^ie,T^-tT,t) , (6.5) 

9^{f,T,t) = g^{Q ^-i'y,T ^ iT,t) . ( 6 . 6 ) 

Generalizing this relation to gauge-invariant dipole-dipole scattering described in 
terms of Wegner-Wilson loops, the IR divergence known from the case of Wegner- 
Wilson lines vanishes and no hnite IR regulator T is necessary. Thus, the Minkowskian 

^In fact, exactly light-like trajectories (7 ^ 00 ) are considered in most applications of the func¬ 
tional integral approach to high-energy collisions [8,11-13,32-34,40-47]. A detailed investigation 
of the more general case of hnite rapidity 7 can be found in [47]. 
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(b) 



Figure 10: High-energy dipole-dipole scattering in the eikonal approximation represented 
by Wegner-Wilson loops in the fundamental representation of SU{Nc)'. (a) space-time and 
(b) transverse arrangement of the Wegner-Wilson loops. The shaded areas represent the 
strings extending from the quark to the antiquark path in each color dipole. The thin tube 
allows us to compare the field strengths in surface Si with the field strengths in surface 
5 *2. The impact parameter b± connects the centers of light-cone momenta of the dipoles. 
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S'-matrix element (6.2), given by the expectation values of two Wegner-Wilson loops, 
forming an hyperbolic angle 7 in Minkowski space-time, can be computed from the 
Euclidean “^'-matrix element” 


{Q,b±,zi,fi±,Z 2 ,r 2 ±) = lim 


{WrACl]Wr,[C2])E 
{WrACl])E{WrAC2])E 


(6.7) 


given by the expectation values of two Wegner-Wilson loops, forming an angle 0 G 
[0, tt] in Euclidean space-time, via an analytic continuation in the angular variable 


S^r2i'y ~ ^^[s/m^],b±,zi,fi±,Z2,f2±) = -ij,b±, zi,ri±, Z2,r2±) , ( 6 . 8 ) 


where E indicates Euclidean space-time and the QCD VEV’s {■ ■ .)e represent Eu¬ 
clidean functional integrals that are equivalent to the ones denoted by (.. .)g in the 
preceding sections, i.e. in which the functional integration over the fermion helds 
has already been carried out. 


The angle 0 is best illustrated in the relation of the Euclidean ^'-matrix ele¬ 
ment (6.7) to the van der Waals potential between two static dipoles, Wir 2(0 = 
0 , b, zi,fi, Z 2 , rg), discussed at the end of this section. 


b±, zi, ri_L, 2 : 2 , fgu) = lim exp 

T —>-oo 


T Wir2(0, b±, ^ 1 , fix, 2 : 2 , r2±) 


(6.9) 


Figure 11 shows the loop-loop geometry necessary to compute 7 ^ 0) ‘) and 

how it is obtained by generalizing the geometry relevant for the computation of 
the potential between two static dipoles (0 = 0): While the potential between two 
static dipoles is computed from two loops along parallel “temporal” unit vectors, 
ti = t 2 = (0,0, 0,1), the Euclidean ^-matrix element (6.7) involves the tilting of 
one of the two loops, e.g. the tilting of ti by the angle 0 toward the X 3 axis, 
ti = (0, 0, — sin 0, cos 0). The “temporal” unit vectors f are also discussed in 
Appendix B together with another illustration of the tilting angle 0. 


Since the Euclidean S'-matrix element (6.7) involves only conhgurations of Wegner- 
Wilson loops in Euclidean space-time and Euclidean functional integrals, it can be 
computed directly on a Euclidean lattice. With (6.7) evaluated numerically for 
many different values of 0 G [0,7r], one needs to hnd the function that describes 
the angular dependence obtained. If this function is analytic in 0, the analytic 
continuation 0 ——iy leads immediately to the desired Minkowskian S-matrix el¬ 
ement (6.2). An obvious difficulty in this proposal is the breaking of rotational 
invariance by the lattice. Moreover, Erst attempts in the direction described have 
shown that the signal size for (6.7) decreases signihcantly with increasing 0 so that 
it is already covered for small values of 0 by the statistical fluctuations [39]. At 
present, it is not clear how to overcome these technical difficulties but the stakes 
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Figure 11: The loop-loop geometry necessary to compute £'^^2 (© / 0)' ‘ ‘) illustrated as a 
generalization of the geometry relevant for the computation of the van der Waals potential 
between two static dipoles (0 = 0). While the potential between two static dipoles is 
computed from two loops along parallel “temporal” unit vectors, ti = t 2 = (0,0,0,1), 
the Euclidean F-matrix element (6.7) involves the tilting of one of the two loops, e.g. the 
tilting of ti by the angle 0 toward the axis, ti = (0,0, — sin 0, cos 0). 


are high: Once precise results are available, the analytic continuation ( 6 . 8 ) could 
allow us to access hadronic high-energy reactions directly in lattice QCD, i.e. within 
a non-perturbative description of QCD from hrst principles. 

More generally, the presented gauge-invariant analytic continuation ( 6 . 8 ) makes 
any approach limited to a Euclidean formulation of the theory applicable for inves¬ 
tigations of high-energy reactions. Indeed, Meggiolaro’s approach has already been 
used to access high-energy scattering from the supergravity side of the AdS/CFT 
correspondence [36], which requires a positive dehnite metric in the dehnition of 
the minimal surface [37], and to examine the effect of instantons on high-energy 
scattering [38]. 

Let us now perform the analytic continuation explicitly in our Euclidean model. 
For the scattering of two color dipoles in the fundamental representation of SU(Nc), 
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the Euclidean ^'-matrix element becomes with the VEV’s (2.14) and (2.35) 


:= m ( 0 , &l, n±, ^ 2 , r 2 ±) 


= lim 

T^oo 


N, + l 

2Nr 


■ exp 


N,-l 


2Nr 


-XS1S2 


+ 


iVc-1 

2Nr 


exp 


iV, + l 

2Nr 


-XSIS 2 


( 6 . 10 ) 


where XSiSj, dehned in (2.27), decomposes into a perturbative (P) and non-perturbative 
(AP) component according to our decomposition of the gluon held strength corre¬ 
lator (2.42): 


,, _ , ,P I ...AP* ^,P _i_ P/AP nc I , .VP c\ 

\S 1 S 2 ^SiS2 ' ^SiS2 ^SiS2 ' \^5'i5'2 ' ^5'iS'2/ 

In the limit Ti = T 2 = T ^ 00 and for 0 G [0, tt], the components read 

xis2 = ^otQx^, xZr = ^otQx^^^\ XS; = cot 0 

with 

(In, - f 2 g\) + g'^D'P (|fig- -r 2 q\) 

- {\riq-f2q\)- g'^Df'^ {\riq - r2q\) 


( 6 . 11 ) 

( 6 . 12 ) 


X^ = 


X 


NP nc 


X 


NP c 


7r^G2(l - p) 
3(iV2 - 1 ) 

7r^G2K 


dT^ (In? - ^25!) + {\nq -r 2 q\) 

- {\nq - r2q\) - {\riq " ^2^-1) 


(6.13) 


(6.14) 


3(iV2 - 1) 


(n-P 2 ) / dvi / (iU 2 (|n <7 + Pin±-P 2 q-P 2 P 2 ±|) (6.15) 


as derived explicitly in Appendix C with the minimal surfaces illustrated in Fig. 11. 
In Eq. (6.13) the shorthand notation g‘^D'p^\\Z 2 _\) = g‘^(\Zi_\) D'p^\\Z 2 _\) is used 
with g‘^{\Zj_\) again understood as the running coupling (2.47). The transverse 
Euclidean correlation functions 

d^K 




(27r) 


D.(K‘)S(K,)S(K^) 


(6.16) 


are obtained from the (massive) gluon propagator (2.43) and the exponential corre¬ 
lation function (2.50) 


Df\zl) = ^Ko{mG\Z^ 
D[^^\zl) = Tra^fd + d'^'' 


a 

Z< 


£){ 2 )(^ 2 ) ^ 27ra^^l-|-^jexp^ 


exp( — 

IZ, 


a 


(6.17) 

(6.18) 
(6.19) 
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With the full 0 dependence exposed in (6.12), the analytic continuation ( 6 . 8 ) reads 

X5i52 = cot0x 0 ^-fy cot(-i 7 )x s^oo ix (6.20) 
-^ ^ 

and leads to the desired Minkowskian S'-matrix element for elastic dipole-dipole 
{DD) scattering in the high-energy limit in which the dipoles move on the light 
cone: 


lim S')^^(s, 6 i,^i,ri_L, 2 ; 2 ,r 2 _L) := \im Z 2 ,f 2 ±) 


= S'f^(cot 0 ^ i, b±, zi,fi±, Z2, f2±) 


= lim 

T—>oo 


iVc+l 

2 W 


exp 


-t- 


.iVe-1 

2 W 


-X 


+ 


N,-l 

2Nr 


exp 


.N, + l 
* 2 W 


-X 


, ( 6 . 21 ) 


where x = X^ + X'^^ + x^^ with (6.13), (6.14), and (6.15). 

It is striking that exactly the same result was obtained® in [ 8 ] with the alternative 
analytic continuation introduced for applications of the SVM to high-energy reac¬ 
tions [11-13]. In this complementary approach the gauge-invariant bilocal gluon 
held strength correlator is analytically continued from Euclidean to Minkowskian 
space-time by the substitution —> —g^p and the analytic continuation of the 

Euclidean correlation functions to real time D^{Z^) D^^z"^). In the subsequent 

steps, one hnds {W[C])m = 1 due to the light-likeness of the loops and that the lon¬ 
gitudinal correlations can be integrated out (Wr^[Ci]Wr 2 [C 2 ])m = f{s,b±, • ■ ■). One 
is left with exactly the Euclidean correlations in transverse space that have been 
obtained above. This conhrms the analytic continuation used in the earlier LLCM 
investigations in Minkowski space-time [8,32-34] and in all earlier SVM applications 
to high-energy scattering [11-13,40-48]. 

In the limit of small y functions, \x^\ -C 1 and 1, Eq. (6.21) reduces to 


lim &i, zi, ri_L, 2 : 2 , r 2 ±) ~ 


1 + 


N!-l 

8iV2 


X 


1 + 


C'2(N,) 2 
4W ^ 


( 6 . 22 ) 


The perturbative correlations, (x^)^; describe the well-known two-gluon exchange 
contribution [93,94] to dipole-dipole scattering, which is, of course, an important 
successful cross-check of the presented Euclidean approach to high-energy scatter¬ 
ing. The non-perturbative correlations, (x^^)^; describe the corresponding non- 
perturbative two-point interactions that contain contributions of the conhning QCD 

®To see this identity, recall that (1T[C']) = 1 for light-like loops and consider in [8] the re¬ 
sult (2.30) for the loop-loop correlation function (2.3) together with the y function (2.40) and 
its components given in (2.49), (2.54), and (2.57) with the transverse Minkowskian correlation 
functions (2.50), (2.55), and (2.58). Note that all these equation numbers refer to Ref. [8]. 
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string to dip ole-dip ole scattering. We analyzed these string contributions system¬ 
atically as manifestations of confinement in high-energy scattering reactions in our 
previous work [33]. 

From the small-y limit, one sees that the full S'-matrix element (6.21) describes 
multiple gluonic interactions. Indeed, the higher order terms in the expansion of 
the exponential functions ensure the fundamental S'-matrix unitarity condition in 
impact parameter space as discussed in [8,45]. 

Concerning the energy dependence, the S'-matrix element (6.21) leads to energy- 
independent cross sections in contradiction to the experimental observation. Al¬ 
though disappointing from the phenomenological point of view, this is not surpris¬ 
ing since our approach does not describe the explicit gluon radiation needed for 
a non-trivial energy dependence. However, based on the S'-matrix element (6.21), 
a phenomenological energy dependence can be constructed that allows a unihed 
description of high-energy hadron-hadron, photon-hadron, and photon-photon reac¬ 
tions and an investigation of saturation effects in hadronic cross sections manifesting 
the S'-matrix unitarity [8,32,34]. This, of course, can only be an intermediate step. 
For a more fundamental understanding of hadronic high-energy reactions in our 
model, gluon radiation and quantum evolution have to be implemented explicitly. 


Although the scattering of two color dipoles in the fundamental representation 
of SU{Nc) is, of course, the most relevant case, we can derive immediately also 
the Minkowskian S'-matrix element for the scattering of a fundamental {D) and an 
adjoint dipole (“glueball” GB) in the Euclidean LLCM. Using (2.40) and proceeding 
otherwise as above, we hnd in the high-energy limit 

lim S'^ge(s, 6 , 2 : 1 , fi, 2 : 2 , f 2 ) := lim S'j!f j^ 2 _i( 0 , 6 , 2 : 1 , n, 2 : 2 , f 2 ) (6.23) 


lim 

T—>00 


iV2- 

^ ’ C 


exp 


■.N, ■ 

l - Y 

. 2 


+ 


A^c+2 
2(iV,+ l) 


exp 


1 ■ 
i-y 

2 


+ 


iVc-2 

2(iV,-l) 


exp 


■ 1 ■ 
i-y 
. 2 


where x = -f- x^^*^ with (6.13), (6.14), and (6.15). 


Finally, we would like to comment on the van der Waals interaction of two color 
dipoles, which is, as already mentioned, related to the Euclidean S'-matrix element 
in the limiting case of 0 = 0 as can be seen from (6.9); The QCD van der Waals 
potential between two static dipoles can be expressed in terms of Wegner-Wilson 
loops [95,96] 


Vr^r2{^ = 0,6, 2:1 = 1/2, fi, 2:2 = 1/2, f2) = 


1 {WrAClWrAC^]) 
{WrACl]){Wr,[C2]) 


(6.24) 


In this limit (0 = 0) intermediate octet states and their limited lifetime become im¬ 
portant as is well known from perturbative computations of the QCD van der Waals 
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potential between two static color dipoles [95-97]: Working with static dipoles, i.e. 
infinitely heavy color sources, there is an energy degeneracy between the interme¬ 
diate octet states and the initial (final) singlet states that leads for perturbative 
two-gluon exchange to a linear divergence in T as T —cx). This IR divergence can 
be lifted by manually introducing an energy gap between the singlet ground state 
and the excited octet state and thus a limit on the lifetime of the intermediate octet 
state [95-97]. 

In the perturbative limit of ^ 0 and T large but finite, i.e. 1) fhe 

perturbative component of our model describes the two-gluon exchange contribution 
to the van der Waals potential which is plagued by this IR divergence due to the 
static limit. In the more general case of finite and T — oo, which is applicable 
also for the non-perturbative component of our model, one cannot use the small-y 
limit and multiple gluonic interactions become important. Here, our perturbative 
component describes multiple gluon exchanges that reduce to an effective one-gluon 
exchange contribution to the van der Waals potential whose interaction range (oc 
I/ttig) contradicts the common expectations. Indeed, it is also in contradiction to 
our results for the glueball mass Mgb which determines the interaction range (oc 
1 /Mgb) between two color dipoles for large dipole separations. As already mentioned 
in Sec. 2.3, we find for the perturbative component, = 2mG, i.e. half of the 

interaction range of one-gluon exchange, by computing the exponential decay of the 
correlation of two small quadratic loops for large Euclidean times r ^ oo, 


Mgb 


y 1. (nf(0)f“/(r)> 

T "{P“»(0)>{P,f(r)> ■ 


(6.25) 


Note that we find for the non-perturbative component, Mq^ = 2/a, which is smaller 
than Mq^ = 2mG with the LLCM parameters and thus governs the long range 
correlations in the LLCM. 

Thus, for a meaningful investigation of the QCD van der Waals forces within our 
model, one has to go beyond the static limit in order to describe the limited lifetime 
of the intermediate octet states appropriately. This we postpone for future work 
since the focus in this work is on high-energy scattering where the gluons are always 
exchanged within a short time interval due to the light-likeness of the scattered 
particles and the finite correlation lengths. Nevertheless, going beyond the static 
limit in the dipole-dipole potential means going beyond the eikonal approximation 
in high-energy scattering and it is, of course, of utmost importance to see how such 
generalizations alter our results. 


45 



7 Conclusion 


We have introduced the Euclidean version of the loop-loop correlation model [8] in 
which the QCD vacuum is described by perturbative gluon exchange and the non- 
perturbative stochastic vacuum model [21], This combination leads to a static quark- 
antiquark potential with color Coulomb behavior for small and conhning linear rise 
for large source separations in good agreement with lattice QCD results. We have 
computed in the LLCM the vacuum expectation value of one Wegner-Wilson loop, 
{Wr[C]), and the correlation of two Wegner-Wilson loops, (Wr^[Ci]Wr 2 [C 2 ]), for 
arbitrary loop geometries and general representations r(j) of SU{Nc). Specifying 
the loop geometries, these results allow us to compute the static quark-antiquark 
potential, the glueball mass, the chromo-held distributions of static color dipoles, 
the QCD van der Waals potential between two static color dipoles, and the S'-matrix 
element for high-energy dipole-dipole scattering. 

We have applied the LLCM to compute the potential and the chromo-electric 
helds of a static color dipole in the fundamental and adjoint representation of 
SU{Nc). The formation of a conhning color hux tube is described by the non- 
perturbative SVM correlations [25] and the color Coulomb held is obtained from 
perturbative gluon exchange. We have found Casimir scaling for both the perturba¬ 
tive and non-perturbative contributions to the chromo-electric helds in agreement 
with recent lattice data [20]. String breaking is described neither for sources in 
the fundamental representation nor for sources in the adjoint representation, which 
indicates that in our approach not only dynamical fermions (quenched approxi¬ 
mation) are missing but also some gluon dynamics. Transverse and longitudinal 
energy density prohles have been provided. For small dipoles, R = 0.1 fm, pertur¬ 
bative physics dominates and non-perturbative correlations are negligible. For large 
dipoles, i? > 1 fm, the non-perturbative conhning string dominates the chromo¬ 
electric helds between the color sources. The transition from perturbative to string 
behavior takes place at source separations of about 0.5 fm in agreement with the 
recent results of Liischer and Weisz [26]. The root mean squared radius R^s of the 
conhning string and the energy density in the center of a fundamental SU{3) dipole 
= 0) are governed completely by non-perturbative physics for large R and 
saturate as R increases at R^°° ~ 0.55fm and = 0) 1 GeV/fm^. 

We have presented the low-energy theorems [28-30], known in lattice QCD as 
Michael sum rules [27], in their complete form in continuum theory taking into ac¬ 
count the important contributions found in [28,31] that are missing in the original 
formulation [27]. We have used the complete theorems to compare the energy and 
action stored in the conhning string with the conhning part of the static quark- 
antiquark potential. The comparison shows consistency of the model results and 
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indicates that the non-perturbative SVM component is working at the renormaliza¬ 
tion scale at which j3{g)/g = —2 and = 0.81. Earlier SVM investigations along 
these lines have found a different value of = 0.57 with the pyramid mantle choice 
for the surface [25,31] but were incomplete since only the contribution from the 
traceless part of the energy-momentum tensor has been considered in the energy 
sum. 

A Euclidean approach to high-energy dipole-dipole scattering has been estab¬ 
lished by generalizing Meggiolaro’s analytic continuation [35] from parton-parton 
scattering to gauge-invariant dipole-dipole scattering. The generalized analytic con¬ 
tinuation allows us to derive A-matrix elements for high-energy reactions from config¬ 
urations of Wegner-Wilson loops in Euclidean space-time with Euclidean functional 
integrals. It thus shows how one can access high-energy reactions directly in lattice 
QCD. First attempts in this direction have already been carried out but only very 
few signals could be extracted, while most of the data were dominated by noise [39]. 
We have applied this approach to compute in the Euclidean LLCM the scattering 
of dipoles at high-energy. The result derived in the Minkowskian version of the 
LLCM [8] has been exactly recovered including the well-known two-gluon exchange 
contribution to dipole-dipole scattering [93,94]. This confirms the analytic contin¬ 
uation of the gluon field strength correlator used in all earlier applications of the 
SVM to high-energy scattering [11-13,40-48]. 

The A-matrix element obtained in our approach has already been used to investi¬ 
gate manifestations of the confining QCD string in high-energy reactions of photons 
and hadrons [33] but leads to energy-independent cross sections in contradiction to 
the experimental observation [8]. The missing energy dependence is disappointing 
but not surprising since our approach does not describe explicit gluon radiation 
needed for a non-trivial energy dependence. In our previous work we have intro¬ 
duced a phenomenological energy dependence into the S'-matrix element that allows 
a unified description of hadron-hadron, photon-hadron, and photon-photon reactions 
and respects the S'-matrix unitarity condition in impact parameter space [8,32,34]. 
However, for a more fundamental understanding of hadronic high-energy reactions 
in our model, one faces the highly ambitious task to implement gluon radiation and 
quantum evolution explicitly. 

More generally, the presented Euclidean approach to high-energy scattering makes 
any method limited to a Euclidean formulation of the theory applicable for investi¬ 
gations of high-energy reactions. Here encouraging new results have been obtained 
with instantons [38] and within the AdS/CFT correspondence [36] and it will be 
interesting to see precise results from the lattice. A promising complementary Eu¬ 
clidean approach has been proposed in [98] where the structure functions of deep 
inelastic scattering at small Bjorken x are related to an effective Euclidean field 
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theory. Here one hopes that the limit a; —0 corresponds to critical behavior in the 
effective theory. The aim is again to provide a framework in which structure func¬ 
tions can be calculated from first principles using genuine non-perturbative methods 
such as lattice computations. In another recent attempt, the energy dependence of 
the proton structure function has been related successfully to critical properties of 
an effective near light-cone Hamiltonian in a non-perturbative lattice approach [99]. 
It will be interesting to see further developments along these lines aiming at an 
understanding of hadronic high-energy scattering from the QCD Lagrangian. 
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A The non-Abelian Stokes Theorem 


In this appendix we review briefly the derivation of the non-Abelian Stokes theo¬ 
rem [16] and explain the emerging surface ordering. We follow the lucid presentation 
given in [10]. 

Let us consider a surface S in Euclidean space-time with boundary C = dS 
and the QGD Schwinger string X;G) defined according to (2.5) that starts 

at some point X on the boundary and evolves along the path C back to the point 
X as illustrated in Fig. 12. We now explain how the non-Abelian line integral over 
C associated with the QGD Schwinger string is transformed into the non-Abelian 
surface integral over S which involves the surface ordering Vs- 

First, we choose an arbitrary reference point O on the surface S and draw a fan- 
type net on S' as a spider could do® (cf. Fig. 12). This net is spanned over S and is 
given by the following curve: Cxo running from X to O, followed by running 

®Note that a real spider draws its net in a sequence different from the one described. The final 
result however is very similar. 
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Figure 12: A surface S with boundary C = dS in Euclidean space-time, the reference 

point O on S, and the fan-type net with center O spanned over S. 

from O to Zi, where the path around the inhnitesimal small square i.e. 

the plaquette at Zi, is attached before it goes back to O along Czio and so on. 
The net is completed with Cxq that runs from O to X. Apart from the initial and 
hnal elements of the net, Cxo and Cxq, we have many plaquettes with “handles” 
connecting them to O. With the following basic properties of the QCD Schwinger 
string, 

^r{0,X-,Cxo)MO,X-,Cxo)-^ = k, (A.l) 

$,(0,X;C'xo)"' = ^liO,X-,Cxo) = MX:0-,Cox = Cxh) : (A.2) 

^riZ,X-,Cxz)MX,0-,Cox) = MZ:0-,Cox + Cxz) , (A.3) 

one sees immediately that the QCD Schwinger string along the net spanned on S is 
equivalent to the QCD Schwinger string <l>r(X, X;C) along the path C: 

<hr(X, X;C) = $,,(0, X; Cxo) ■ (pi’oduct of QCD Schwinger strings 

for the plaquettes with handles) ■ <Fr(0,X; Cxo)~^ ■ (A.4) 

Next, we consider the contribution of a single plaquette with handle. The QCD 
Schwinger string for one plaquette, say the one at Zn singled out in Fig. 12, reads [10] 

(plaquette at Z^) = 1 - ig^Aa^^{Zn) Q%{Zn) + ■ ■ ■ . (A.5) 
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where Aa^y{Zn) denotes the surface element at the point Zn- Taking into account 
the handles, the contribution of this plaquette to (A.4) becomes 

$^(0, Zn] (plaquette at Z„)<h^(0, Z^] Cz„o)~^ 

= 11 - ig^Aa,,{Zn) g%{0, Z^] Cz^o) K+--- , (A.6) 

where (A.2) and the parallel transported gluon field strength as defined in (2.4) have 
been used. 


Finally, inserting (A.6) into (A.4) for all n, i.e. summing up the contributions of 
all plaquettes with handles, while respecting the ordering, one obtains in the limit 
of an inhnitesimally hue net 


^r{X,X]C)=^r{0,X]Cxo) 


(A.7) 


■ Vs exp 



da^,{Z)G<^ 


flV 


{0,Z]Czo)K 


^r{0,X]Cxo)-^ . 


Here Vs denotes the ordering on the whole surface S as implied by the net shown in 
Fig. 12. Taking the trace in (A.7) and exploiting its cyclic property leads ultimately 
to the non-Abelian version of Stokes theorem 


Tr^ <l>r{X,X]C) 


Tr^ Vs exp 



da^n{Z)G%{0,Z]Czo)t. 


(A.8) 


B Loop and Minimal Surface Parametrizations 


A rectangular loop Q with “spatial” extension Ri and “temporal” extension 2Tj 
placed in four-dimensional Euclidean space, as shown in Fig. 13, has the following 
parameter representation: 


with 





a = ^ 

cfucf 

U Cf U Cf 


(B.l) 

= ] 

[xfi 

[Ui) 

= Ao, 

; - (1 - Zi] 

)ri + Uiti, M, 


\ (B.2) 

Cf = < 

[xf- 


= Ao, 

; - (1 - Zi] 

)ri + Viri + Ti ti, 

n, e [0,1]] 

[ (B.3) 

C? = \ 

'xf( 

[Ui) 

= Ao, 

[ + ZiTiV', 

Ui ti , tii 

e [T„-T,]] 

\ (B.4) 

cf = -1 

[xF- 


= Ao, 

; - (1 - a! 

)ri + Viri- Ti ti, 

n, e [1,0]] 

\ (B.5) 
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Figure 13: (a) Spatial arrangement of a color dipole and (b) its world-line in Euclidean 
“longitudinal” space given by the rectangular loop Ci that defines the minimal surface Si 
with dSi = Ci- The minimal surface is represented by the shaded area. In our model, it 
is interpreted as the world-sheet of the QCD string that confines the quark and antiquark 
in the dipole. 
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where 


Ti 


^ Ri sin6'j cos0i ^ 
Ri sin^j sin0j 
Ri cos 9i cos0j 
y Ri cos 9i sin 0j 


and ti 


( ° ^ 
0 

— sin &i 
\ cos Oi y 


(B.6) 


The “center” of the loop Q is given hj Xq^. The parameters Zi^ Ri, 9i, (pi, and 
0j are dehned in Fig. 13 which illustrates (a) the spatial arrangement of a color 
dipole and (b) its world-line Ci in Euclidean “longitudinal” space. The tilting angle 
0j 7^ 0 is the central quantity in the analytic continuation presented in Sec. 6. 
Moreover, 0i = 7r/2 together with 02 = 0 allows us to compute conveniently the 
chromo-magnetic held distributions in Appendix C. 


The minimal surface Si is the planar surface bounded by the loop Q = dSi given 
in (B.l). It can be parametrized as follows 

Si = ^^Xi{ui,Vi) = Aoi - (1 - Zi)ri + Viri + uRi, Ui G [-Tj,Tj], n* G [0,1]| (B.7) 

with Ti and given in (B.6). The corresponding inhnitesimal surface element reads 


d<7^i/ (A^j) 


dXj^ dXiu 
dui dvi 


dX,^ dXj, 
dvi dui 


dui dvi 






dUi dVi . 


(B.8) 


C X Computations with Minimal Surfaces 


The quantities considered in the main text are computed from the VEV of one loop 
{W[C]) and the loop-loop correlation function (yV[Ci\W[C 2 \)■ Using the Gaussian 
approximation in the gluon held strengths, both are expressed in terms of XSiSj 
functions (2.15) and (2.27) as shown in Secs. 2.1 and 2.2. These y functions are 
central quantities since here the ansatz of the gauge-invariant bilocal gluon held 
strength correlator and the surface choice enter the model. In this appendix, these 
functions are computed explicitly for minimal surfaces (B.7) and the ansatz 

given in (2.42), (2.44), and (2.48). Note that the contributions from the inhnitesi- 
mally thin tube—which allows us to compare the held strengths in surface Si with 
the held strength in surface S 2 —cancel mutually. 

Depending on the geometries and the relative arrangement of the loops, the y 
functions determine the physical quantities investigated within the LLCM such as 
the static qq potential (3.2), the chromo-held distributions of a color dipole (4.9), 
and the S'-matrix element for elastic dipole-dipole scattering (6.10). 
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We compute the three components xfiSj’ Xsisl separately for general 

loop arrangements from which the considered quantities are obtained as special 
cases. Without loss of generality, the center of the loop C 2 is placed at the origin 
of the coordinate system, Xq 2 = (0,0, 0,0). Moreover, C 2 is kept untilted, ©2 = 0, 
and 0 := ©1 is used to simplify notation. We limit our general computation to 
loops with ri ^2 = (hi, 2 ±, 0 , 0 ) = 6 * 1^2 = 7 r /2 and transverse “impact parameters” 
b = Xqi — Xq 2 = dfoi = ( 6 i,& 2 , 0 , 0 ) = ( 6 j_, 0 , 0 ) which allows us to compute all of 
the considered quantities. 


X^^S2 Computation 


Starting with the dehnition 


NPc 
A Si 52 


TT 


4 Js; 

Ti‘^G2H 


da,,{X,) I dap^iX2) = W - X 2 ) 


'S 2 


(C.l) 


10/'Ar2 / d(Jpcr{X2) {6fj_p6i/cr dpadup) D(^Z ) , 

~ -*-1 JSi JS 2 

one exploits the anti-symmetry of the surface elements, da^^ = —da^p, and applies 
the surface parametrization (B.7) with the corresponding surface elements (B. 8 ) to 
obtain 


xZs2 = ^ose 


with 


2 11 

.NPc _ „„„ Q f . in ■ r 2 ) [ dvi [ dv 2 [ dui [ du 2 D{Z‘^) 

~ 1 ) Jo Jo J-Ti J-T 2 

/ b±- { 1 - Zi) fi_L Vi fi_L + (1 - Z2) r21. - V2 r2_L \ 

— Ui sin 0 

\ Ml COS 0 — M2 / 


(C.2) 


Z = Xi - X 2 = 


, (C.3) 


where the identities ti ■ r 2 = ri ■ t 2 = 0 and ti ■ t 2 = cos 0, evident from (B. 6 ) with 
the mentioned specihcation of the loop geometries, have been used. In the limit 
T 2 —cx), the U 2 integration can be performed 


cT2 


lim 

T2^oo 


du2D{Z^) = 


-T 2 

d^K 

d^K 


d^K 


('T 2 


D{K^) lim / du 2 e 


iKZ 


T2^oo 


I-T 2 


D{K‘^)2'k 5{Ki) exp 


d- iK^ui sin 0 -|- iK^ui cos 0 


^(3) (^2 


exp 


iK±Z^ + iK^ui sin 0 




(C.4) 
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which leads to 


lim Xs^sl = cos 0 


n'^G2K 


{ri±-r2±) [ dvi f dv2 f duiD^^\Z‘^) . (C.5) 

Jo Jo J-Ti 


3(iV2-l) 

Taking in addition the limit Ti — oo, the ui integration can be performed as well 

[ = i 

-Ti 1 limTi^oo2Ti 


lim 

Ti—>-oo 


27r6{K3 sin 0 ) = 


I sin©I 


for sin 0 7 ^ 0 , 
for sin 0 = 0 . 


With Ti = T 2 = T/2—i^cx), one obtains for sin 0 7 ^ 0 


jvpc COS 0 7 r^G' 2 K 


lim - 

iim XS1S2 - 


T—^00 

and for sin 0 = 0 

NPc 


sin 0 | 3(iV2 — 1 ) 
ti‘^G2H 


{G±-r2±) / dvi / dv2D^‘^\Zl) 


lim X 5 TS 2 = lim T cos 0 —(fi± ■ / dvi / dv 2 D^^\Z^) 

T^oo T^oo 3(iV^^-ij Jo Jo 


(C. 6 ) 


(C.7) 


(C. 8 ) 


Evidently, (C.7) is the result given in ( 6 . 12 ) and (6.15) which describes the conhning 
contribution to the dipole-dipole scattering matrix element S'pp. 


From (C.8), one obtains the conhning contribution to the static color dipole 
potential for Si = S 2 = S which implies Ti = T 2 = T/2, 0 = 0, zi = 2 : 2 , ri = r 2 = r, 
and ri_i_ ■ r 2 j_ = so that 

lim Xss^ = lim T / dvi f dv 2 D^^\Z^ = {vi - V 2 f R^) 

T^oo T^oo 3(7V2 _ ij 

^ ^ [ dp{l- p)D^^\p^R^) , (C.9) 

T^oo - Ij Jo 

which leads directly to (3.7). 


From (C.5) the conhning contribution to the chromo-held distributions AG^^(X) 
can be computed conveniently. Equation (C.5) reads for Si = S'p, Ti = Rp/2 and 
Ri = Rp, and S 2 = Sw, T 2 = T/2 and R 2 = R 


tt'^G k G G r^pG ^ 

lim = COS0 (Tip ■f2p) / dvi / dv2 / duiD^^\Z^) (C.IO) 

3(iV^ - Ij Jo Jo J-Rp/2 

with 


Z = Xi - X 2 = 


6 j_ - (1 - Zi) fi± + Vi fi_L + (1 - Z 2 ) f2± - V 2 r2± 
ui sin 0 


. (C.ll) 
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The confining non-perturbative contribution to the chromo-magnetic fields vanishes 
as it is obtained for plaquettes with 0 = 7r/2. The corresponding contribution to the 
chromo-electric fields can be computed with 0 = 0 as follows: Due to Ri = Rp ^ 0, 
the Ml and Vi integrations in (C.IO) can be performed with the mean value theorem. 
Keeping only terms up to 0{Rl), the confining non-perturbative contribution to the 
chromo-field distributions is obtained as given in (4.18). 


Computation 


We start again with the definition 


TT 


xZsT ■■= J^dap^X^) J^dap^{X,)F/:j;:^{Z = X,-X,) 

^ 12(4 - ij 


(C.12) 


1 

X- 

2 


^ {ZJpp-ZpSp,) + ^{ZpS, 


dZ, 


ua ' ^crdup) 


Di{Z^ 


and use the anti-symmetry of both surface elements to obtain 
^2 

XsiS2 


NPnc _ ^G2{1 k) f f dap,{X2)^Z^6ppD,{Z^) 

JS2 

r 

da,UXi)^ [da^{X,_)j^S„D[(Z^ 


with 


D[{Z^) = 


6(iV2-l) Js, 

7r‘^G2{l-K) 

WFXJs; 

7r^G2(l-fi;) 

■ 3(iV2 - 1) Jsl 

"d^K 


lvJS2 


dx. 


2(7 




d 

dK^ 


Di{K^) . 


(C.13) 

(C.14) 

(C.15) 

(C.16) 


As evident from (C.15), Stokes theorem can be used to transform each of the surface 
integrals in XsiS^'^ ^ integral: 


NPnc 

AS 1 S 2 


7r2G2(l n) f^ 


3(iVc^-l) Js, 

7r^G2(l-fi;) 


dZ, 


dZp{X2)6ppD[{Z^ 


iyJC2 


6(iv2 -1) jsr^^ 

7r^G2(l-K) 


da^,(Xi) (b dZp{X2) 6pp Z, D,{Z^) 


C 2 


3(^2 _ 1 ) jJ^ZpiX,) jJiZp{X2) 6pp D[{Z^^ 


(C.17) 

(C.18) 

(C.19) 
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With the line parametrizations of Ci and C 2 given in (B.l) and the specihcation of 
the loop geometries mentioned at the beginning of this appendix, (C.19) becomes 


NPnc _ _ ^^^ 2(1 1 ^) 

AS 1 S 2 - _ X) 


(C.20) 


r rn 

I‘T2 


< cos 0 / dui 

/ dU2 

D[(ZlA-D[{Zlc)-D'AZhA)+D\{Zlc) 

y J-Ti j 

'-72 



+ (nx • f2i)fdv, fdv2\D[(ZlB)-D[{Zlo)-D[{ZlB)+D[(Zlr,) 

Jo Jo 


where the following shorthand notation is used: 


Zxy := Wf - X 2 ^ with G and X 2 ^ G C. 




X 


yX 




■yy 


(C.21) 


In the limit R 12 ^ 2 — 00 , the term proportional to ri_L-r 2 x on the RHS of (C.20) 

can be neglected and thus (C.20) reduces to 


r-Ti 


r*T2 


lim Xa^o'^ = ~cos 0—lim / dui lim / du 2 (C.22) 
l\ZZ 3(iV2-l) ^ ’ 


X 


D((Z1^) - D[{Zic) - + D[{Zl.c) 


Here, the integrations over ui and U 2 can be performed analytically proceeding 
analogously to (C.4) and (C. 6 ). With Ti = T 2 = T/2 ^ 00 , one obtains for 
sin 0 7 ^ 0 


lim Xsi 


NP nc 
S2 


COS0 7r^G2(l — k) 

sin 0 | 3{N^ — 1) 


X 


D'J\z'UJ - D'nzlcJ - D'P'(Zl„^) + B'X'(Zjcx) 


y(2) 


'( 2 ), 


h(2)/ 


(C.23) 


and for sin 0 = 0 


lim XsfsT ^ ~ lim T cos 0 ^ ^ 

3(iV2-l) 


(C.24) 


X 


d'JHzIa) - Dnz'ic) - D'r(ziA) + D'r(zic) 


V ( 3 ) / '72 


’ ( 3 ) / ^2 


■)' ( 3 ) ( ry‘^ 


With the identities 


ZaAL — ^Iq — fgg , ZacI. — ^ig — r2q , ZqAX — ” fgg , Zcc± — ” fgg , 

(C.25) 
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one sees immediately that (C.23) is the result given in (6.12) and (6.14) that de¬ 
scribes the non-conhning non-perturbative contribution to the dipole-dipole scatter¬ 
ing matrix element S^d- 


From (C.24), one obtains the non-conhning contribution to the static potential 
for Si = 82 = S, i.e. Ti = T 2 = T/2, 0 = 0, ri = r 2 = r, 


hm = - hm T f 

T^oo 3(iV2-l) 


X 


Df’iMA) - Df\zlc) - Df\zl^) + 

which contributes to the self-energy of the color sources with 

7r^G2(l - «) 


lim = 

XS 5 self - 3 (iV 2 - 1 ) 


lim T 


Df\zl^) + D’r{Z'l,,.) 


d3). 


= - li,„ T 

T^oo 3(iV2-l) ^ ^ 

and to the potential energy between the color sources with 


lim xTs;:t 

1 -^00 


= hm T 


7r^G'2(l - k ) 
6 (iV| - 1) 
7r^G2(l - «) 


£>;'='(ZX) + D'8(Zl,) 

T-;.- 3(7V-1) 

The last gives the non-conhning contribution to the static potential (3.6). 


= hm T 


(C.26) 


(C.27) 


(C.28) 


The non-conhning non-perturbative contribution to the chromo-electric helds 
[AG2^(X) with af3 = 14, = 4i\ can be computed most conveniently from (C.18) 
with zero plaquette tilting angle 0 = 0. The corresponding contribution to the 
chromo-magnetic helds [AG^^(X) with af3 = ij = ji] is obtained for plaquette 
tilting angle 0 = 7r/2 and thus vanishes which can be seen most directly from the 
surface integrals (C.13). Now, we set 0 = 0 to compute the contribution to the 
chromo-electric helds: Using the surface Si = Sp and loop U 2 = dSw parametriza- 
tions, (B.7) and (B.l), with our specihcation of the loop geometries, one obtains 
from (C.18) 


nc 

X.Sp S\y 


to 

to 

1 = 

1 

l'Rp /2 

/ dui 

J—Rp /2 J 

[ dVi 

0 

3(A2 - 1) 




< / dU2 

(0± ■ 

)r>i(Zj^)-(ru.Zici)Ci(Z?c) 

[ J-T/2 

*- 

-* 


(C.29) 


- {ri± ■ r 2 i_) 



(r[_i_ ■ ZiB±) Di{ZIb) — (hi± ■ Zipp) Di{Z\b) 
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with Ti = Rp/2, Ri = Rp, T 2 = T/2, R 2 = R, and the shorthand notation 

Zix ;= Xi - Xf with Xi e = Sp and Xf e = dS^ . (C.30) 

In the limit R T ^ cx), the term proportional to ■ r 2 ± on the RHS of (C.29) 
can be neglected, 


NPnc 


limxspSw 


7t‘^G2{1-k) 

3(iV| - 1) 


cRpR 


dui 


'-Rp/2 


dvi lim 

T —>00 



X 


(ri_L ■ ZiAA_) — (ri_L ■ Zic±) Dii^ic) 


(C.31) 


With an inhnitesimal plaquette used to measure the chromo-electric held, Ri = 
Rp —> 0, the mean value theorem can be used to perform the ui and vi integrations 
in (C.31). Keeping only terms up to 0{R^), this leads directly to the non-conhning 
non-perturbative contribution to the chromo-held distributions AG^^(X) as given 
in (4.16) and (4.17). 


Computation 


Comparing the dehnition of the perturbative component 


"ZS\S 2 


TT 


g 


/ da,.{X,) / dapM2) F^^pAZ = - X2) 

ISl JS 2 


= — / d(jp,AXi) I dap AX 2 ) 

A 


(C.32) 
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1 

X- 

2 


^ {ZJpp-Zp5pA + -^{Zp5, 


dZ, 


ua ~ Z/jdyA 


Dp{Z‘^ 


with that of the non-conhning non-perturbative component Xs^AA (C. 12 ), 

one hnds an identical structure. Thus, accounting for the diherent prefactors and 
the diherent correlation function, the results for x^iS 2 read oh directly from 

the results for XsAA above. 

With Ti = T 2 = T/2 ^ 00 and our specihcation of the loop geometries, one 
obtains the result for sin 0 7 ^ 0 from (C.23) 


lim xfiS2 = - 


cos 0 
I sin 01 


(C.33) 


X 


r>p"’(zLi) - D'p’izlcA.) - r>7'(zSAi) + D'p"’{zhci) 


y(2)/ 


y(2)/ 


h(2)/ 
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(C.34) 


and the result for sin 0 = 0 from (C.24) 


lim XsrS 2 = - 1™ T cos Qg^ 

1 —^oo 1 —^OO 


X 


- -Dp'"'(Ziio) - + Bp'''(Zjc) 


y(3)/ 


V(3)/ 


V(3)/ 


where Zxy is dehned in (C.21) and Zxy± is given explicitly in (C.25). Evidently, 
(C.33) is the hnal result given in (6.12) and (6.13) which describes the perturbative 
contribution the dipole-dipole scattering matrix element Sdd- 


The perturbative contribution to the static potential is obtained from the ex¬ 
pression corresponding to (C.26), 


lim Yqq = — lim T 

T-,oo T^oo 


X 


DfHzl.) - DTizlc) - c;!“'(zS.4) + D'^‘'{Zlc) 


JL) 


y(3)/ 


which contributes to the self-energy of the color sources with 


lim Xs5self = 

1 —>oo 


lim T 

T^oo 


D 


'(3), 


p (ZL) + r>p”(zJc) 


= - lim T 2 D'^^\z 


T—>oo 


AA) 


and to the potential energy between the color sources with 


(C.35) 


(C.36) 


lim Xsspot = - lim T g 

1 —^oo 1 —>oo 


Df\z'ic) + DY'(ZiA) 


'(3), 


Mm T 2g‘ Df\Z\c) ■ 


T—>oo 


(C.37) 


The last gives the perturbative contribution to the static potential (3.4). 


The perturbative contribution to the chromo-magnetic helds [AG^^(X) with 
a(5 = ij = ji\ vanishes while the one to the chromo-electric helds [AG^^(X) with 
aP = i4 = 4i], for which a plaquette with 0 = 0 is needed, is obtained from the 
expression corresponding to (C.31), 


\im XspSw 



r‘Rpl 2 


dui 


<-Rp/2 


dvi lim 

T—>oo 



X 


if L± ■ Ziay) DpiZ‘lf) — (ri_i_ ■ Zic±) DpiZ^f 


(C.38) 


with Zix as dehned in (C.30). To perform the Ui and Vi integrations in (C.38), 
again the mean value theorem can be used since the plaquette has inhnitesimally 
small extensions, Ri = Rp ^ 0. Keeping only terms up to OiRp), this leads directly 
to the perturbative contribution to the chromo-held distribution AG^^(A) as given 
in (4.13) and (4.14). 
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